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Abstract 

We present a new method for computing the Konishi anomalous dimension in A/" = 4 SYM at weak 
coupling. It does not rely on the conventional Feynman diagram technique and is not restricted to 
the planar limit. It is based on the OPE analysis of the four-point correlation function of stress- 
tensor multiplets, which has been recently constructed up to six loops. The Konishi operator 
gives the leading contribution to the singlet SU{4) channel of this OPE. Its anomalous dimension 
is the coefficient of the leading single logarithmic singularity of the logarithm of the correlation 
function in the double short- distance limit, in which the operator positions coincide pairwise. 
We regularize the logarithm of the correlation function in this singular limit by a version of 
dimensional regularization. At any loop level, the resulting singularity is a simple pole whose 
residue is determined by a finite two-point integral with one loop less. This drastically simplifies 
the five-loop calculation of the Konishi anomalous dimension by reducing it to a set of known 
four-loop two-point integrals and two unknown integrals which we evaluate analytically. We 
obtain an analytic result at five loops in the planar limit and observe perfect agreement with the 
prediction based on integrability in AdS/CFT. 
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1 Introduction 



It has been realized recently that the four-point correlation function of the so-called stress-tensor 
multiplets in A/" = 4 super- Yang-Mills theory (SYM) has a new symmetry [1]. In combination 
with A/" = 4 superconformal symmetry, it imposes strong constraints on the integrand of the loop 
correction to the correlation function and leads to an iterative structure at weak coupling, at any 
loop order and for a gauge group of arbitrary rank [2]. The correlation function of four stress- 
tensor multiplets plays a special role in A/" = 4 SYM theory. In virtue of the operator product 
expansion (OPE), its asymptotic behaviour at short distances contains information about the 
anomalous dimensions of a large variety of Wilson operators and the corresponding structure 
constants of the OPE. Moreover, if considered in the planar limit and restricted to the light cone, 
it is dual to the four-particle scattering amplitudes [3l S]. 

For more than ten years, this correlation function was not known beyond two loops. The 
main difficulty in going to higher loops is due to the factorially increasing number of contributing 
Feynman diagrams. In general, each individual diagram respects neither gauge invariance nor 
conformal symmetry, but the symmetries are restored in the sum of all diagrams. This calls 
for developing a new approach that makes full use of A/" = 4 superconformal symmetry and of 
the specific symmetry of the stress-tensor multiplet mentioned above. Such an approach has 
been proposed in two recent papers [H |2] , where a new construction of the four-point correlation 
function was carried out in A/" = 4 SYM theory for the gauge group SU{Nc) with arbitrary Nc- In 
the planar limit, for iVc — )■ oo and with the 't Hooft coupling a = g'^Nc/ (47r^) fixed, the integrand 
of the four-point correlation function was determined up to six loops and the non-planar 0{1/N^) 
correction was identified at four loops (up to four arbitrary rational constants). □ 

In this paper, we apply the results of Refs. [H [2] to perform the OPE analysis of the four- 
point correlation function of the stress-tensor multiplets. As the main result of our analysis, we 
present a new method for computing the Konishi anomalous dimension in A/" = 4 SYM theory for 
arbitrary gauge group SU{Nc). The Konishi operator is the simplest unprotected gauge invariant 
Wilson operator in A/" = 4 SYM, whose scaling dimension receives anomalous contribution to all 
loops. In the OPE context, the distinguishing feature of the Konishi operator is that it controls 
the leading asymptotic behaviour of the four-point correlation function at loop level in the short 
distance limit. In this manner, we obtain an analytic result for the Konishi anomalous dimension 
at five loops in planar A/" = 4 SYM theory and observe perfect agreement with the prediction 
based on integrability in AdS/CFT El E]- 

The properties of the Konishi operator have been studied extensively after the discovery of the 
so-called Konishi anomaly [9] in supersymmetric gauge theories. The interest in the subject was 
renewed in the context of the AdS/CFT correspondence. As was observed in [10], the Konishi 
supermultiplet in A/" = 4 SYM theory is a long (or unprotected) multiplet that corresponds 
to the first string level in the spectrum of type IIB excitations in an AdSs x background. 
Recently, the Konishi anomalous dimension 7k: (a) again attracted a lot of attention after the 
discovery of integrability in the planar limit on both sides of the AdS/CFT correspondence (for 
a review see [H]). At strong coupling, the first few terms of the expansion of 7Ac(fl) in powers of 
^-1/4 jj^ |;]^g planar limit were obtained from the semiclassical quantization of short strings on an 
AdSs X background [121 [131 El]- At weak coupling, the values of 7j<;(o) at four and five loops 
in planar Af = 4 SYM were predicted in Refs. [S] and [B] , respectively, from the integrable string 

^It is worth mentioning that our approach is not limited to six loops. Extending it to higher orders is just a 
question of computer power. 
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sigma model by evaluating finite size effects using Liischer's formulasH The four-loop prediction 
was later confirmed by direct perturbative calculations using JV = 1 Feynman super-graphs 
[T^ [TB] and ordinary Feynman diagrams [T7] . Until now, no five-loop test of the integrability 
prediction had been performed, and it is not clear whether traditional techniques would allow 
one to reach such a high perturbative level. A numerical prediction for 7x;(a) at intermediate 
coupling, interpolating between the strong and weak coupling results, was made in [18] from the 
solution of the F— system of integral non-linear equations and more recently in fL9\ from the 
TBA equations. 

As was already mentioned, the Konishi operator provides the leading contribution to the 
asymptotic behaviour of the four-point correlation function at short distances G(l,2,3,4) ~ 
(^12)''^'^''"^^^ as Xi — )■ X2- At weak coupling, this asymptotic behaviour implies that perturbative 
corrections to the correlation function at i loops are given by a sum of logarithmic singularities 
(In 0:^2)'^ with powers k < i. The coefficients of the higher powers of logarithms (for k > 1) are 
expressed in terms of the anomalous dimensions at lower loops. It is only the single logarithm 
(with k = 1) that carries information about the anomalous dimension at i loops. This fact 
complicates the evaluation of the anomalous dimension. It is more advantageous to consider 
instead the logarithm of the correlation function lnG(l, 2, 3, 4), whose asymptotic behaviour at 
short distances involves a single logarithmic singularity to all loops. We can further simplify 
the analysis by considering the double short- distance limit Xi — )■ X2 and X3 — )■ X4, in which case 
lnG(l,2,3,4) ~ (7;c(a)/2)[lnx I2 + In 0:34] . To determine the anomalous dimension in this way, 
we need an efficient way of computing the perturbative corrections to lnG(l, 2, 3, 4). 

Applying the results of Ref. 0, we can express lnG(l,2,3,4) at £ loops as a Euclidean 
integral, whose integrand is a conformally covariant function of the four external points and the 
£ integration points. In the short-distance limit Xi — )■ X2, the integral develops a single logarithmic 
singularity ~ lnx^2 when all £ integration points approach the external point Xi simultaneously. 
Similarly, for X3 — )■ X4 the logarithmic singularity ~ lnx|4 originates from integration in the 
vicinity of the point X3. It is clear that these singularities are of ultraviolet (UV) origin with 
the small distances XI2 and x'^^^ playing the role of a UV cut-off. To extract the coefficient of 
the logarithmic singularity of the integral, which defines the anomalous dimension of the Konishi 
operator, we can simplify the calculation by taking X12 = X34 = inside the integral and by 
introducing the most convenient regularization scheme for the resulting UV divergences. This 
is done by changing the integration measure from D = AtoD = 4 — 2e dimensions. Thus, 
we transform the expected single logarithmic singularity of lnG(l,2,3,4) in the double short- 
distance limit into a simple pole 1/e. Our final simplification comes from the observation that, 
for xi = X2 and X3 = 0:4, the £— loop residue at this simple pole is in fact given by an (£ — 1)— loop 
finite two-point integral of the propagator type. We can then apply an array of well-known and 
very efficient methods for computing such integrals. 

In summary, we have reduced the problem of computing the Konishi anomalous dimension 
at £ loops to the problem of evaluating a finite two-point integral at (£ — I) loops^ This allowed 

^See also [3 [8] for an alternative approach using the mirror thermodynamic Bethe Ansatz. 
•^This feature is typical for usual renormalization group calculations in momentum space |20) . 
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us to obtain the following result up to five loops @: 



lic{o,) = 3a — 3a^ + — a' 
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where the non-planar four-loop correction is predicted up to an arbitrary rational constant, r. 
We use (n to denote values of the Riemann zeta-function ({n) at integer points. The obtained 
expression for ■jjc{a) is in agreement with the existing perturbative four-loop results of [151 CHI 
[T7t and with the five-loop prediction of |51 El El E] based on integrability in AdS/CFT. As 
a byproduct of the OPE analysis of the four-point correlation function, we also investigated the 
spectrum of anomalous dimensions of the twist-two operators with non-vanishing spin at three 
loops and found agreement with the values conjectured in [23]. 

The paper is organized as follows. In Section 2 we define the four-point correlation function of 
stress-tensor multiplets in A/" = 4 SYM and use the OPE to relate 7Ac(a) to the leading asymptotic 
behaviour of the correlation function in the short distance limit. In Section 3 we formulate our 
method for computing the Konishi anomalous dimension and illustrate it by evaluating 7/c(a) up 
to two loops. In Section 4 we extend our analysis to four loops. Applying the results of Refs. [2] 
for the four-loop correlation function and making use of well-known techniques for evaluating 
Feynman integrals, we express the four-loop correction to '~fic{0') as a linear combination of six 
master two-point three-loop integrals (five in the planar sector and only one in the non-planar 
sector). If rewritten in the dual momentum representation, the latter coincide with some known 
finite three-loop integrals of the propagator type. In Section 5 we evaluate the five-loop correction 
to jjcio-) in the planar limit. We show that it is given by a linear combination of 22 master scalar 
four-loop integrals. Among them 20 integrals correspond to planar graphs and coincide, in 
the dual momentum representation, with known finite four-loop propagator integrals [27]. The 
remaining two non-planar integrals are evaluated in Appendix [Bl Section 6 contains concluding 
remarks. In Appendix A, we describe the method of IR rearrangement in the configuration space 
that we employ in our calculation of 7ac(o). In Appendix C, we perform the OPE analysis of 
the four-point correlation function and extract the values of three-loop anomalous dimensions of 
twist-two operators with Lorentz spin zero, two and four. 

2 Four-point correlation function 
2.1 Expression for the integrand 

In this paper we study the OPE of the stress-tensor multiplet in A/" = 4 SYM. This is the simplest 
example of a half-BPS operator, whose superconformal primary state has the form 



''The one- loop value of 7)c was found for the first time in [21]. At two loops, it was first extracted from the OPE 
of two stress-tensor multiplets in A/" = 4 SYM, as part of the investigation of the two-loop four-point correlation 
functions of half-BPS operators [22l[23]. The three- loop value, together with the anomalous dimensions of all 
twist- two Wilson operators in A/" = 4 SYM was originally conjectured in [21]. This three- loop prediction for 7^; 
was then confirmed for the first time in [25 by a direct perturbative calculation. 



0^^, = tr($^<l>^)- Vtr ($^$^) . 
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It is built from the six real scalars $^ (with / = 1,...,6 being an 5*0(6) index) in the ad- 
joint representation of the gauge group SU (Nc) and belongs to the representation 20' of the R 
symmetry group 5*0(6) ~ 5f/(4). To keep track of the 50(6) tensor structure of the OPE, it 
proves convenient to introduce auxiliary 50(6) harmonic variables Yj, defined as a (complex) 
null vector, Y"^ = YjYj = 0, and project the indices of O^"^ as follows: 

0{x,y) = YiYjOiiix) = YiYj tr {x)^\x)) , (2.2) 

where y denotes the dependence on the y— variables. 

An important property of the operator (12. ip is that its scaling dimension is protected from 
perturbative corrections. The same is true for the two- and three-point correlation functions of 
the operator 0{xi,yi). The four-point correlation function is the first to receive perturbative 
corrections: 

oo 

G4= (0(xi,yi)0(x2, y2)O(x3,l/3)0(x4, 1/4)) =5Za^ of (1,2,3,4), (2.3) 

e=o 

where the expansion on the right-hand side runs in powers of the 't Hooft coupling a = g'^N^/ (47r^) 

and of denotes the perturbative correction at i loops. Notice that here we do not assume the 

(I) 

planar limit and allow G4 to have a non-trivial dependence on Nc. It is this four-point function 
that will serve as the starting point of our OPE analysis. 

At tree level, Cf^ reduces to a product of free scalar propagators and the corresponding 
expression can be found in Ref. [I]. At loop level, the superconformal symmetry of the Af = 4 
SYM theory restricts Of to have the following factorized form [I] '■ 

of (1,2,3,4) = ^^^^ X /?(1,2,3,4) x ^^(a;,), (for £ > 1) , (2.4) 

where F^^\xi) is a function of Xi only (with i = 1, 2, 3, 4) to be specified below and i?(l, 2, 3, 4) is 
a universal, £— independent rational function of the space-time, Xi, and harmonic, Yj, coordinates 
at the four external points 1, 2, 3, 4, whose explicit expression can be found in Ref. [I]. 

So from (12.31) and (12. 4p we find that the loop corrections to the four-point correlation function 
are determined by a single function F^^\xi). As was shown in Refs. [H |2], this function has a 
number of remarkable properties in A^ = 4 SYM theory. Namely, it can be represented in the 
form of an loop Euclidean integral, 

,-yj2 ^2 ^2 ^2 ^2 ^2 p 

d\,...d^x,^,f'\x,,...,x,^,), (2.5) 

where the integrand f^^^ depends on the four external coordinates xi,...,a;4 and the ^ addi- 
tional (internal) coordinates xs, . . . ,^4+^ giving the positions of the Lagrangian insertions. The 
integrand Z*^^-* can be written in the form 



/W(a;i, . . . ,0:4+.) = - v-'---'-:-^ . (2.6) 



pW(a;i,...,a;4+£) 



Hi 



<i<j<4:+l ^ij 



Here the denominator contains the product of all distances between the (4 + £) points and P*-^^ is 
a homogeneous polynomial in xjj of degree {i — !){£ + 4)/2. Most importantly, this polynomial 
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is symmetric under the exchange of any pair of points Xi and xj (both external and internal). As 
we have shown in Refs. |2j, this property alone combined with the correct asymptotic behaviour 
of the correlation function in the short- distance and the light-cone limits, allows us to completely 
determine F^^\xi) up to six loops in the planar sector. In the non-planar sector, the same analysis 
leads to an expression depending on a few constants only. For example, at one and two loops, 
we have 

p(l) _ 1 p(2) _ -'- 222_222, (O 7\ 

~ ; ~ 43 2-^ ^0-10-2 ^0-3 0-4^0-50-6 ~ ^12^34^56 + • • • ) l^-' j 

o-gSe 

where in the second relation the sum runs over all 5*6— permutations of the indices 1, . . . , 6. 
Similar expressions at higher loops can be found in Ref. [2]. 



2.2 Operator product expansion 

As was mentioned in the previous subsection, the four-point correlation function (12. 5p has a 
particular asymptotic behaviour at short distances, dictated by the operator product expansion 
(OPE). 

For the scalar operators (12.21) the OPE takes the following form 



Oix^,y^)Oix2,y2) = c^ ^^^'J'^ X + cM ^ ^2) + co^^^^O'^^,{x2) + ... (2.8) 

X]_2 \^12) ' ^12 



where we only displayed the contribution of operators with naive scaling dimension up to two. 
Here the most singular l/xf2 contribution comes from the identity operator X while the first 
subleading 0(l/a;^2) contribution originates from two operators: the half-BPS operator (12. ip 
and the Konishi operator defined as 

/C = tr (<l>^$^) . (2.9) 

Since the operators X and C20' are protected, the constants cj and co do not depend on the 
coupling constant and keep their tree-level values, cj = (A^"^— l)/(327r'^) and co = l/(27r^). For the 
Konishi operator, both the coefficient cic{a) and its scaling dimension A^(a) receive perturbative 
corrections to all loops. In what follows we will mostly concentrate on the anomalous dimension 
of the Konishi operator 

00 

A;c = 2 + 'jjcia) = 2 + ^ a^^'^ . (2.10) 

e=i 

In the singular limit Xi — )■ X2 we can apply the OPE expansion (12. 8 p to find the asymptotic 
behaviour of the correlation function (12. 4p at short distances (in Euclidean kinematics). It 
receives contributions from all operators on the right-hand side of (12. 8p . A crucial advantage of 
the Konishi operator is that it has the minimal possible scaling dimension among all unprotected 
operators. To separate the contribution of the Konishi operator it is useful to consider the 



^The Konishi operator is the simplest of an infinite series of t-wist-t-wo operators contained in the OPE 
In Appendix C we extract from the OPE the three-loop anomalous dimensions of the t"wist-two operators "with 
Lorentz spin 2 and 4. 
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double short- distance limit xi — )■ X2, X3 — )■ x^. Taking into account the relation fl2.8p we obtain 
the asymptotic behaviour of the four-point correlation function in this limit as 



^3 {N^-lfyf.yl, 
' 4(47r2)4 xf^xi^ (47r2)4 



?/?2l/34(l/?3?/|4 + l/lVi3) 



2 2- 

rf-'-^ /y^ 

1 nJ^oA -] 



+ ^ Jlllr (4(a)«-(")/^ - 1) 



12-^ 34 -^13 

■ \ ... , (2.11: 



O Ai2"''34"''13 

where m is a conformal cross-ratio defined in fl2.13p below, yf- = [Yi ■ Yj) denotes the scalar 
product of harmonic variables and the dots denote subleading terms. 

Comparing the OPE prediction (12.111) with the general expression for the correlation function, 
Eqs. (12. 3p and (12 ■4p . we obtain the following relation for the functions F^^\xi) for X2 — )■ xi and 

X4 — 7- 

^a^FW(x.) -V (4(a)u^'=(")/^ - 1) X [1 + 0{u) + 0(1 - v)] . (2.12) 

Here u and v are the two conformally invariant cross-ratios made of the four points x^, 

99 99 
^ = ^12^, ^ = ^i^, (2.13) 

13"''24 "''13"''24 

so that u — )■ 0, f — )■ 1 in the double short-distance limit X2 — )■ , 3:4 — )■ 0:3. For our purposes it 
is convenient to introduce the notation for the function xf^F^^\xi) in this limit, 

xtsF^'\x.) ^3 Fix.), (2.14) 
and to rewrite the OPE limit (I2.12p as 

ln(^l + 6Y,a'F^'\xi)^ ^ In u + In {cl{a)) + 0{u) + 0{{1 - v)) . (2.15) 

Let us now expand both sides of the relations (I2.12p and (I2.15P in the powers of the coupling 
a and compare their short-distance asymptotics. We find from (I2.12p that F^^\xi) ~ (Inw)^ as 
M — > 0. In particular, from (I2.12p we have to two-loop order 

F« = l7i:^)ln« + V) + ..., 

F(^) = 1(7«)^ (Inu)^ + (^7i?) + l^'M') + ^«(^) + . . . , (2.16) 

where the constants 7^"* and a^^-* define the perturbative corrections to the anomalous dimension 
7a;('j^) = Yle>i^^'yK^ coefficients (c^(a))2 = 1 -|- 3 a^a^^^. In a similar manner, 

from (I2.15P we obtain that the particular combination of the functions F'^^\xi) with 1 < < ^, 
arising from the expansion of the logarithm on the left-hand side of (I2.15p . scales as Inw. For 
instance, at two-loop order we have from (I2.15P 

F(^) - 3 (F«)^ = ^7^^ \nu + \a^'^ - ^(a«)^ + . . . . (2.17) 
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Comparing this relation with fl2.16p . we observe that the two-loop correction to the anomalous 
dimension appears in (12.161) in the subleading term, while in (I2.16P it defines the leading 
singular behaviour. As we show in the next section, this property can be used to drastically 
simplify the calculation of the Konishi anomalous dimension 7a: (a). 



3 Method for computing the Konishi anomalous 
dimension 

Here we present our method for computing the Konishi anomalous dimension at higher loops. It 
takes full advantage of the known properties of the correlation function explained in the previous 
section. In this section we illustrate the key features of the method with the simplest examples 
of one and two loops. 

Before we do this, we would like to recall the standard approach for extracting the Konishi 
anomalous dimension from the asymptotic logarithmic behaviour of the four-point correlation 
function (see, e.g., Refs. [23 123] )• With the help of the relations (12. 5 p - (12.71) we obtain the 
following expressions for the correlation function to two loops: 

F«=^?(l,2,3,4), 

= h{l, 2; 3, 4) + h{3, 4; 1, 2) + h{2, 3; 1, 4) + h{l, 4; 2, 3) 

+ h{l, 3; 2, 4) + /i(2, 4; 1, 3) + ^ {xj^xl^ + + ^lA) 2, 3, 4)]^ . (3.1) 

Here the notation was introduced for the one- and two-loop conformal Euclidean integrals 

1 [ (i^X5 



ly*^ nn" /y*^ nn" 

"''15"''25"''35"''45 



h{\ 2- 3 41 = I ^^^5 c?^a^6 

z,, o, '±; 2^2 / (^2 ^2 ™2 \™2 {^2 ™2 ™2 \ ' y^'^l 

with the remaining /i— integrals obtained by permuting the indices. 

The explicit expressions for these integrals as functions of the conformal ratios (I2.13P are 
known [30], but what we need here is just their asymptotic behaviour for x\ — )■ x^ and X3 — )■ X4, 
or equivalently m — )• and f — )• 1. Replacing the integrals in (13. ip by their asymptotic expansions, 
we easily obtain the following result for the one- and two-loop correlation functions in the singular 
short- distance limit: 

F« = lln«-1 + ..., 

F(2) = l(lnn)2-ln« + ^C3 + ^ + ---, (3-3) 



leading to 



F(2)-3(F«)^ = -ilnn + ^C3 + l + ... • (3.4) 



These relations are in perfect agreement with the OPE prediction (I2.16P and (I2.17p . They allow 
us to reproduce the well-known result for the two- loop Konishi anomalous dimension, Eq. (II. ip . 
and the two-loop normalization coefficients, a*-^-* = —1 and a^^-* = 3^3/2 + 7/2. 
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3.1 One loop 

Let us now return to the one-loop expression F'^^\ Eq. (13. 3p . and understand the origin of the 
singularity F^^'^ ~ Inu at short distances. It is easy to see from (13.21) that for xi — )■ X2 and 
— )■ X4 the integral 5^(1, 2, 3, 4) develops a logarithmic divergence coming from the two distinct 
integration regions x^ — )■ Xi and — )■ X3, 



F« L / L / (3 5) 



(i*^X5 1 f (i^Xs 

Here we have restricted the integration to two balls of radius 6, centred at the points Xi ~ X2 
and X3 X4. Choosing xf^ ^ 5^ ^ xfg, 3^34 allows us to replace the other two propagator factors 
in the first integral X35 ~ X45 by ^^3, and similarly for the second. This simplification of the 
integrand by replacing it with its asymptotic expression in the relevant integration region will be 
very helpful at higher loops. Going to radial coordinates we find 

(3.6) 

where the short distances x\2,x1^ — )■ serve as UV cut-offs and the dots denote terms finite in 
the limit x\2-,x\^ — > 0. It is easy to see that this relation is in agreement with the first relation 
in ( 1331) . 

Let us now examine what happens if we interchange the integration in (13. 5 p with taking the 
limit Xi — )■ X2, — )■ X4. In this limit, the first integral on the right-hand side of (13. 5 p reduces 
to d'^x^/x\i and it diverges as 0:5 xi. This is not surprising since x\2 plays the role of 

a short-distance cut-off in the first integral in (13. 5p and (13. 6p . Therefore, in order to define the 
integral for = we have to introduce a different short- distance regulator. The simplest way 
to do this is to modify the integration measure in (I33|) as follows % 

rf^xs ^ /i^' (f-'^'x5 , (with e < 0) , (3.7) 

without changing the form of the integrand. In this way, we find from (13.51) 

9-'-^-^ I ^-^/ (3.8) 



where we introduced the subscript in F^^^ to indicate that it is defined in the regularization 
scheme (13. 7p at x\2 = = 0. Going to spherical coordinates, d'^~^''x = S^r^~'^^dr with = 
27r^~^/r(2 — e), we find that the two integrals in (13. 8p produce equal contributions, leading to 

Fi^) = + 0(e°) = 1 + 1 ln(/.V^^) + . . . . (3.9) 

Comparing the right-hand sides of the relations (13. 6p and (13. 9p . we observe that they coincide 
(up to the 0(l/e) term) upon the identification — ^ /^^ and — )■ /i^. In other words, for 



^We would like to emphasize that this regularization is different from the conventional dimensional regulariza- 
tion in coordinate space in the sense that we modify the integration measure only and use the scalar propagators 
instead of l/(x^)^~'^. This explains why e should be kept negative. 
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^12 — ^34 — 0' within the regularization scheme (13. 7p . the dimensionful parameter yu^ plays the 
role of the UV cut-off. This property allows us to relate the coefficient in front of Inu in the 
asymptotic behaviour of F^'^^ at small u, Eq. (13. 3p . to the residue of Fe^^ at the simple pole 
1/e. Moreover, it follows from (I2.16P that this coefficient coincides with the one-loop Konishi 
anomalous dimension 7^^-', leading to 

in agreement with (II. ip . 

This suggests a new method for computing the Konishi anomalous dimension: Instead of 
evaluating the finite four-dimensional integrals in (13. ip and finding their asymptotics at u — )■ 
0,v ^ 1 afterwards, we can first evaluate the integrand at xi = X2 and 0:3 = X4, thus making 
the integrals divergent, then introduce the regularization (13. 7p and, finally, identify the terms 
singular for e — )■ 0. 

We would like to emphasize that this method captures correctly only the terms divergent 
for M — > 0, f — )■ 1 but not the finite ones. To see this, let us apply the above procedure to the 
one-loop expression F^^^ 

,2e r ^4-2e^ ^4 



47r2 J x\^xl^ 



In comparison with ( 13. Sp here we did not restrict the integration region over X5. This is not really 
necessary, since the integral converges at large X5. To perform the integration, it is convenient 
to switch to the dual momenta k = X15 and p = X13. Then, the integral in (13. lip takes the form 
of the standard one-loop "bubble" momentum integral of propagator type: 

M-.^-i^!!/,^, ,3.12) 



47r2 J k\p-kY 
leading to (with Jj? = n^/{e^^ vr)) 

= M(^) = {xlJu'r + ^ + 0(e^)) . (3.13) 

Here the pole in e comes from the integration over small momenta A; — )■ and (j9 — fc) — )■ and, 
therefore, has an IR origin in the dual momentum representation. Comparison of (I3.13P with 
the first relation in (13. 3 p shows that the singular term is reproduced correctly (upon identifying 
^125^34 ~^ /^^ subtracting the pole), while the regular (constant) term is different. 



3.2 Two loops 

Let us now extend the above analysis to two loops. According to (13. 3p . the two-loop correction to 
the correlation function F*^^) has a stronger, (InM)^ singularity for -u — )■ 0. The reason for this is 
that, in the short distance limit xi — )■ X2 and X3 — ?■ X4, the integrals on the right-hand side of (13. ip 
develop overlapping singularities when the integration points X5 and Xq independently approach 
the two external points, e.g. X5 — )• Xi and — )• X3. At the same time, the leading (InM)^ 

^In what follows, for the sake of simplicity we do not distinguish between and ^"^ . 
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singularity is supposed to cancel in the particular combination of one- and two-loop corrections 
(13.41) . which defines the 0(0^) correction to the logarithm of the correlation function on the 
left-hand side of fl^Jf^ . 

To understand the reason for this, we replace F*^^^ and F'^^-' on the left-hand side of fl3.4p by 
their explicit expressions ( 13. ip and, then, simplify the resulting expression by applying the same 
limiting procedure as in (13. lip . Namely, we take the limit xi — )■ X2 and — )■ X4 inside the g— 
and /i— integrals and modify the integration measure as in (13.71) . In this manner, we arrive at 

Fi2) _ 3 (^(i))2 = ( / d^-^-xsc/^-^-xe '^^^^^^/f + ^,^f f .^^^"^^^^ , (3.14) 

v4vrV J {xtAJxiei^b^e) 

where the expression on the right-hand side is manifestly symmetric with respect to the inte- 
gration points, X5 and xq, and it takes into account the contribution from the sum of g^— and 
/i— integrals. It is clear from (I3.14p that the integral diverges logarithmically when X5 and xq 
approach the external points xi and X3. The simplest way to evaluate (I3.14p is by going to the 
dual momenta ki = X15, ^2 = xiq and p = X13, so that 

- 3 (fW)2 = 4M(2) - 2(M«)2 . (3.15) 

Here the one-loop integral M*^^^ was introduced in (I3.12p and M ^^^^ stands for the standard two- 
loop scalar propagator-type integral [3T] 

M(2) = 



fl^'Y f d^-^'kid'-^'k2 



47rV J kfkliki - k2Y{p - kiY{p - 



+ + (3,16) 

Substituting (I3.13p and (I3.16p into (I3.15p . we find that the double pole cancels leading to 

_ 3 (P(l))2 ^ (^2^/^2)-2e (^-1 - ^ + 0(e)) . (3.17) 

We observe that the expansion of this expression around e = produces the logarithmic term 
— (1/2) ln(/i2/x^3), which matches the — (l/4)lnM term on the right-hand side of (13. 4p after the 
identification x\2)X\^ — )■ /i^. As in the one-loop case, the finite terms in the two expressions are 
different. According to (I2.17p . the coefficient in front of Inw is related to the two-loop Konishi 
anomalous dimension, 7^^ /12. Similarly to (I3.10p . this allows us to write 

7?^ = 6^[F(^)-3(Fi^))^]=-3, (3.18) 

in agreement with (II. ip . 

What is the reason why the double pole cancels in (I3.14p ? This happens because the nu- 
merator in (I3.14P has the following characteristic feature: it vanishes for x^ — )■ xi and x^ — )■ X3 
with in general position. As a consequence, the most singular contribution coming from the 
two regions, x^ — )■ Xi^Xq — j- X3 and x^ — )■ x^.Xq — )■ Xi, is suppressed and the integral in (I3.14p 
develops a weaker singularity. It only arises when the two integration points approach one of the 
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external points simultaneously, X5,xe — )■ xi and x^jXq — )■ 0:3. We can make use of this fact to 
further simplify the calculation of the divergent part of the integral (I3.14p . 

Like in the one- loop case fl3.5p . we can single out the divergent contribution to fl3.14p by 
introducing a dimensionful parameter 5^ ^ xf^ and restricting the integration in f l3.14p to a ball 
of radius 6 centred at the points xi or X3. Due to the symmetry of the integral f l3.14p under the 
exchange of xi and X3, the two regions produce the same contribution leading to 



V^TT / Jq^ [Xi^XiQjX^Q 

Here the integration is performed over the region Vl^ defined as xf^jX^^ < S^. Notice that in 
this region we can safely replace X35 and x|g inside the integral with xfg. To simplify f l3.19p . we 
introduce the radial coordinates r| = x^^, = x^i and the angle (p between the two vectors, 
■ Xiq) = r^TQ cos (j). Integration over the angle yields [32] 

r.h.s. of » = 1^ ^J-l = + 0(«»)) , (3.20) 

where the notation was introduced for r< = min(r5,r6) and r> = max(r5,r6). As expected, the 
residue of the pole in fl3.20p is the same as in f l3.17p and, therefore, it leads to the same result 
for the two- loop Konishi anomalous dimension (I3.18p . 

3.3 Loop reduction 

Comparing the integrals in Eqs. (]3.14p and (]3.19p . we notice that the latter contains a smaller 
number of propagators and, therefore, is much easier to analyze. Still, both integrals are two- 
loop ones as they involve integration over two points. As we show in this subsection, there is yet 
another simplification which allows us to effectively eliminate the integration over one point and, 
therefore, reduce the number of loops by one. 

Let us return to the relation fl3.19p and perform the integration over Xg with the point X5 in 
a general position inside the region Qs- It is easy to see that the Xg— integral depends on the two 
scales xf^ and 6^, which play the role of cut-offs at short and large distances, respectively. Notice 
that the integral converges at large Xq and, therefore, it would stay finite if we sent S^/xl^ to 
infinity. This means that, as far as the leading divergence of fl3.19p is concerned, when computing 
the Xg— integral we can neglect its 5^— dependence and extend the integration over xg to the whole 
(4 — 2e)— dimensional (Euclidean) space 

= a{xl,rx[l + 0{xlJ6')]. (3.21) 

Here the second relation follows from dimensional analysis of the integral and is some constant 
regular at e = 0. Substituting this result on the right-hand side of fl3.19p we arrive at the integral 

' ' [l + 0ixy5^)] . (3.22) 



27r2 ./o. ixl,y+^ 
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cx+P-D/2 



Figure 1: Diagrammatic representation of the chain relation. Sohd hne with index a stands for 
l/(a;^)" and black dot denotes the integration point. 

It is easy to see that the ©(xfg/^^) term does not produce a pole in e and, therefore, can be 
safely neglected. Performing the Xs— integration, we finally obtain 

) - 3 ~ + O(e0)^ . (3.23) 

Thus, the residue at the pole and, hence, the two-loop Konishi anomalous dimension, is deter- 
mined by the constant which is given in its turn by the one-loop integral f l3.2ip 

C. = &f = 1 + 0{e) . (3.24) 

Substituting this relation into f l3.23p . we arrive at f l3.20p . 

The key point in the above argument is that the final result in f l3.24p is given by an integral 
which is finite for e — i- 0. Indeed, the potential singularities of this integral could come from 
integration in the vicinity of Xi (at infinity it is obviously finite by power counting). However, 
the numerator of the integrand vanishes for xg — )■ Xi, so the integral is convergent. Notice that 
the relation ( 13.24p can be rewritten in a form resembling the first line in (I3.19P : 



.2 



27r2 



+ 



^16%6 



r d^-^'xe 


f d^-^'XG 






J xIqXIq J 





(3.25) 



In this form, each integral in the square brackets develops a pole 1/e but their sum is finite for 
e — 7- 0. The first integral diverges when xq Xi, the second one when xq — t- oo. Both of them 
can be easily evaluated with the help of the "chain relation" [32] shown diagrammatically in 

Fig.m 

1 f d^xo G{a,(3) 



It is easy to check that the pole 1/e cancels in the sum of two integrals. The third integral in 
( I3.25P is scaleless and, therefore, it vanishes in dimensional regularizationj^ 

3.4 The method 

We are now ready to formulate our method for computing the Konishi anomalous dimension at 
higher loops. It consists of four steps: 



^More precisely, it develops two poles 1/e after integration over xq xi and — ?• oo, one of UV and the 
other of IR origin. Their residues have opposite signs, so they cancel in the sum. 
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Step 1: Expand the logarithm of the correlation function (12.151) in powers of the coupling 
constant and obtain the £— loop correction to the left-hand side of (I2.15P in the form of an 
folded integral over the internal points x^, . . . , x^+f 

In I 1 + 6^a^F(^)(a;i) I = XI / ^^^5 . . .d^x^+ili{xi, . . . ,Xi\x5, . . . ,x^+i) , (3.27) 

\ £>1 J 1>1 

with the integrand symmetric under the 6*4 x 5*^ permutations of the four external points, 
Xi, . . . , X4 and the i integration points X5, . . . , x^+i. In the short distance limit, for xi — X2 and 
X3 — )■ X4, the integral develops a single logarithmic singularity. 

Step 2: Replace the integrand le by its limiting value at xi = X2 and X3 = X4, introduce the 
regularization (13. 7p for the integration measure over the i internal points and, then, restrict the 
integration over X5, . . . , ^4+^ to a ball of radius 6 centred at one of the external points, say, Xi. 

2 X X a^fi^' / d^-'^'x5 . . . d^-'^'xi+i Ii{xi\x5 . . . , x^+e) . (3.28) 

Here we inserted the factor of 2 to take into account the contribution from the integration around 
the point X3, and introduced the notation for 

^lim^ Ie{xu ■ ■ .,3:41x5 . . . ,Z4+£) = Ii{xi\x5 . . .,X4+e) + Oid^/xj^) . (3.29) 

— >x^ 



Step 3: Freeze one of the integration points, say, X5 and perform the integration over the 
remaining {£ — 1) points by extending the integration region to the whole space: 

2 I d'-'^x,...d'-'^x,+Mx,\x,...,x,+,) = ^(x?5)-2-(^-i)% (3.30) 

with the constant C^-i being regular at e = 0. Going back to (13.271) . we perform the remaining 
rcs— integration and obtain: 



In 1 + 6 5^ a^FW(x.) ~ ^'^eJ-^ f 
\ e>i J £>i 



-J2^'C,.iyif^y6T + 0{e'). (3.31) 



The expansion of this relation in the powers of e produces a In /x^ term which should match the 
Inu term on the right-hand side of (12.151) for XI2, — )■ /i^. 

Step 4: We compare the last relation with (I2.15p . identify Inw — )■ 21n(/i^/(5^) and obtain the 
Konishi anomalous dimension as 



7/c(a) 



-^,n(l.6E.'Pi')w) (3.32) 
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leading to 



lda) = -Y,^'C^-i- (3-33) 

i>i 

This relation allows us to express the Konishi anomalous dimension at i loops in terms of the 
constants C^-i which are defined in their turn in terms of scalar integrals at {C, — 1) loops, 
Eq. ^M)- 

As we will show in this paper, going through these steps we will be able to determine the 
Konishi anomalous dimension up to five loops in the planar sector, as well as at four loops in the 
non- planar sector (up to a rational prefactor, see below). 



4 Konishi anomalous dimension at three and four loops 

Let us apply the method described in the previous section to compute the Konishi anomalous 
dimension at three loops and beyond. 



4.1 Preliminaries 



We start with the general expression for the logarithm of the correlation function in the short- 
distance limit on the left-hand side of fl3.27p . 



In [l + eJ^a^Fw] = J]a^/W 



(4.1) 



where /'•^•* 


are given 


by the following expressions up to five loops: 


j(i) 


= 6FW, 




J(2) 


= 6[F(2) 


-3(F«)^], 


J(3) 


= 6[F(3) 


-6F«F(2) + 12(F«)3], 


/(4) 


= 6[F(^) 


_ 6F(i)f(3) - 3(F(2))2 + 36F(2)(f(i))2 - 54(F(i) 


J(5) 


= 6[F(5) 


_ 6F(i)f(^) - 6F(3)f(2) + 36F(3)(F(i))2 



+ 36F«(F(2))2 _ 216F(2)(F«)3 + i|6(P(i))5] 



(4.2) 



We recall that the function F^^^ was defined in fl2.14p as the short-distance limit of the correlation 
function at i loops, Eq. (12. 5p . 

To find we substitute (12. 5 p into (I2.14p and regularize the integration measure according 
to dSH): 



fW(x,1 



(4.3) 



£! ( Arr'^y J {xf^xl^,) . . . {x\ ,^_^_gX^^^j^^) Y[5<i<j<4+e^ij 
Here P^^^ coincides with the polynomial P^^\ Eq. (12. 6p . evaluated at X2 = Xi and X4 = x^ 

PW= lim pW(xi,...,a;4+,). (4.4) 



X4 — >-j^3 
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For instance, for £ = 1 and £ = 2 we find from (12. 7p 



P*- •* — 1 , ^ — 2x^^x^Q + Ax^^x^^x^Q + Ax^^x^qX^^ . (4-5) 

Notice that the polynomial P'-^-* is symmetric with respect to all (4 + i) points while for the 
polynomial P*^^^ this symmetry reduces to 5*2 x 5"^ permutations of two external points Xi, x^ and 
of the i integration points x^, . . . , x^+i. 

Substituting the definition fl4.3p into fl4.2p . we can represent the right-hand side of (14. ip in 
the same form as in fl3.27p . in terms of the £— fold integrals 

= J d'-'^x, . . . d'-'^x,+e X, , (4.6) 

and express the integrands Xi in terms of the polynomials P^^'-* with For i = \ and 

£ = 2 we have 



-7- 3 x\^ 



■^is / -'5.6 ^ 4 I _ J -^Isl-^lS-^Se "T ■^16-^35 ■^13-^56J 

"•^13 — /,_4 ^4 ^4 ^4 ^4 ^2 ' ' ) 



( /\ 'fr 2^2 fy»4 ,-y»4 ryt^ /y»4 l /y»2 1"3 / /] -77*4 <-y»4 /y»4 ,-y»4 /y»4 /-y»2 

l^" J "''15'^35"''16'''36 \ -^56 / ^" •^15-^35-^16-^36"''56 

where Ps^g = -P'-^H^S; ^e)- It is easy to see that the numerator of I2 vanishes for 0:5 ^ and 
2^5 ~^ 2:3, as needed to achieve a lower degree of divergence. 

As was already explained, the integral in (14. 6 p develops a simple pole 0(l/e) from integration 
over all points X5, . . . , ^4+^ in the vicinity of the two external points xi and X3. For — > xi we 
can safely replace — )■ inside Xi without affecting the residue at the pole. In this way, we 
construct the function Xg defined in (I3.29p . At one and two loops we have 

T _ 3 1 _ 3 {xfr^ + xIq — Xgg) , , 

By construction, the functions Xi do not depend on 2:3. Substituting these relations in (I3.30p . we 
find with the help of ^^?2^ 

Co = -3 , Ci = ^C, = ^, (4.9) 

with C, defined in ^2^. 



4.2 Three loops 

According to (12.50 and (12.60 the correlation function at three loops F^^\xi) is determined by the 
5*7— invariant polynomial P^'^\xi, . . . , X7). As was shown in Ref. [2], the form of this polynomial 
can be fixed from the requirement for the correlation function to have the correct asymptotic 
behaviour at short distances. The result is 

P^^-* = — (3^12)^(^34^45^56^67^73) + '^'7 permutations, (4-10) 
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where the sum runs over the permutations of the indices 1, . . . , 7. The exphcit expression for 
P^^) contains 5040 distinct terms. However, as was explained above, for our purposes we only 
need its limit fl4.4p for X2 = Xi and X4 = x^. This brings the number of terms down to 27: 



0(3") o22442i,i422222i/i62222 

V, / >t ly ly rf /-y» /y" _l_ /\/y /y /y y /y ry _l_ /l/y /y y /y /y 

r — OJ'i3J'i6"''l7"f'35"''36 "r ^"''13-^16"''17'''35"''37"''56 ^"''13"''17'''36"''56"''57 

.i3a.^7a.37.x,gQ + permutatious , 



io44 2 2 2io4 2 2 4 2io6 2 2 4 

_i_ / ly T* ni* T* ly —X— / 'y nf y y /y _1_ / nf nf /y nf 



(4.11) 



where the permutations only act on the integration points X5,Xq,Xy. The polynomial P'-^-* 
defined in this way is a completely symmetric function of X5,xq,xi. 

Then, we apply f l4.3p and (14. 2 p to define the integrand in (14. 6 p at three loops: 



X 



13 



(4vr2)'n,=5,6.7^f.4 



p 



5,6,7 



2 2 2 

<-yj ^ /y^ /y ^ 

■^56 -^67 ■^75 



— 6x 



13 I 



X 



^5,6 -^6,7 -^5,7 



72a; 



56 



67 



13 



57 



(4.12) 



where Pjj was defined in (14. 7p and the notation was introduced for P5,6,7 = P*-^-* (^5; ^6, 2:7). 
Replacing the P— polynomials by their explicit expressions we get 



2x' 
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/' /I 'TT- 2 '\ 3 /y^ /y^ <-y>4 /y^ /y^ ly 



67 



06222, 422/22 

-iry ly ly ly 1 ly ly ly 1 ly nr* 

'J"''13"''56"''57'''67 "r •^13-^17"''56l"''37'''56 



10X35X57) 



,2 2/2 2 4,92 2 2 2,2 2 4\,/|2 4 4 2,q 
+^13^56 1^35^36^17 + -^^16^35^37^17 + ^15^16^37^ + ^^16^17^35^36 + pcrm. 



. (4.13) 



We verify that the numerator X3 vanishes when one of the integration points X5, xg, X7 approaches 
the external points Xi or X3. This ensures that the integral in (14. 6 p develops a single pole only. 
It originates from two different integration regions, X5,X6,X7 Xi and X5,X6,X7 — )■ X3, which 
produce however the same contribution in virtue of the symmetry of X3 under the exchange of 
xi and X3. 

We examine X3 for X5,X6,X7 — )■ xi and determine the corresponding function (I3.29P at three 
loops 



X3 



1 



^ /I •77- 2 ^ 3 /-y»4 /y^ /-yi4 ^2 ,-yi2 ^2 



A i 2 . 42, 24io222| 222 

/I /y ly _l_ /y /y _l_ /y ly _1_ / y _l_ T* 'T* 

^•^15-^16 -^IR-^fi? "T 't'lR-^fi? "T ^J'lR't'lfiJ'fiy "T J'lRJ'lfiJ'fi 



'15-^67 



'15-^67 



"15-^16-^67 



a5-^16-^56 



— lOxJgXggXgy + 3x5gX57Xg7 + 5*3 permutatious 



(4.14) 



As before, X3 is a symmetric function of X5,Xg,X7 and it does not depend on X3. At the next 
step, we substitute the function X3 into the left-hand side of (I3.30p and integrate it over Xg,X7 
with X5 fixed 



^^xeci'* ^^X7X3(xi|x5, xe, X7) 



74-2e 



G 



277 



(4.15) 



In close analogy with the situation at two loops (see Eq. (I3.25P and the discussion around it), 
the integral on the left-hand side of ( I4.15P is finite for e — )■ 0, due to the special properties of the 
expression in the square brackets in (14.141) . We would like to emphasize that this only holds for 
the sum of all terms in the square brackets, but not for each individual term. In other words, 
if we split the integral on the left-hand side of (I4.15P into a sum of integrals corresponding to 
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each terms in fl4.14p . then each integral develops a pole 1/e. The poles cancel in the sum of all 
integrals leading to a finite expression for the coefficient C2 in firi5|l . % As an example, let us 
examine the integrals corresponding to the two terms in the second line of (I4.14p : 



4 4 4 2 2 2~ 

/->'■ ^ / >"» ^ /-) ■ ^ (T^*-" (^-^ 



10 [ d 



X 



15 



j4-2e 



Xq 



X' 
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4-2^x7 3 

— H 

4 ^2 ^ ^4 

15 



M-2e 



Xe 



X 



16 



j4-2e 



Xj 
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(4.16) 



where in the second relation we took into account that the integral over xq is scaleless and, 
therefore, it vanishes in dimensional regularization (see footnote E]) • 

Examining the contribution to fl4.15p from the remaining terms on the right-hand side of 
( I4.14p . we find that most of the integrals vanish in the same manner as in fl4.16p . The remaining 
non-zero contribution takes the following form: 



Co 



d^-^'x^d^-'^xj 



4-2e, 



4 2 2 



1 



4 4 2 2 

4 

+ 



2 2 2 2 

/y ^ /y^ /y ^ /y^ 



'Alii 

/y^ /y^ <-yj-" ly-' 

X^gX^^XggXg^ 


' 1 i 1 1 

Xj^gXj^yXggXgy 


+ 


2 2 2 2 

nn'-' ly ^ nn-" ly^ 

X]^gX;^7X5gXg7 


2 


1 




4 


~^ 4 2 2 2 

<-yj^ /-y>-^ <-yj-" ly'-' 


~^ 2 2 2 2 
"''16'^17-^ 56-^57 


+ 


4 2 2 2 
■^16"''56"''57"''67 


4 


1 -^67 


+ 


2X57 


4 2 2 2 2 

/y^ <-yj-" /-yj-" /y-' ly^ 


4 4 2 2 

Xj^gXj^yXggXgy 


2 4 2 2 

<T^ ^ ly^ nn-" ly^ 

X]^gX;^7X5gXg7 



(4.17) 



Here we set x\^ 



1 for simplicity since the dependence of the integral on x\^ is uniquely fixed 
by dimension analysis. It is convenient to represent the 12 integrals in this relation in the form 
of Feynman diagrams as shown in Fig. [2l 



Xi 












Figure 2: Diagrammatic representation of the integrals in f l4.17p . Solid lines without labels 
depict propagators l/a^^j, the label 2 referring to the square of the propagator. Dashed lines 
represent numerator factors xf^, black and white dots represent integration and external points, 
respectively. 



We observe that the first 9 integrals in f l4.17p can be easily evaluated with the help of the 
"chain relation" f l3.26p shown in Fig. [H Applying f l3.26p consecutively, we can express the first 9 



^To be more precise, the distribution (xf^) ^ ^'^ is singular for e 
pole 1/e whose residue is a contact term, see Appendix A. 



0. So, the left-hand side of (|4.15p has a 
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integrals in terms of the G— function and obtain the following relation (with D = A — 2e) 
1 



2G(1, 1)G(2, 2 - D/2) + 2G{1, 1)G(1, 4 - D/2) + 2G(2, 1)G(1, 4 - D/2) 
+2G(1, 1)G(1, 3 - D/2) + (G(2, 1))' + 2G(1, 1)G(2, 1) + (G(l, 1))^ 



+4G(1, 1)G(2, 3 - D/2) + 4G(1, 1)G(1, 3 - D/2) + Jio + In + /; 
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(4.18) 



Here Jio, /n and J12 stand for the last three integrals on the right-hand side of f l4.17p . corre- 
sponding to the last three diagrams in Fig. [2J The simplest way to compute them is to introduce 
the dual momenta ki = xiq, ^2 = Xn and p = X15 and rewrite the above integrals as two- 
loop propagator-type momentum integrals. The latter can be easily computed using MINCER [31], 
yielding 

/io = 4(G(l,l)f [l-l^ + l^' + 0{e') 
/n = (G(l,l))2(-2 + 4e-2e2 + 0(e3)) , 

lu = 2(G(1, 1))2 (1 - 2e + 4e2 + ©(e^)) . (4.19) 

Substituting these relations in fl4.18p we verify that all poles in e cancel on the right-hand side 
of (14.181) . leading to the Konishi anomalous dimension at three loops, Eq. 03.331) . 



in agreement with (11. ip . 



7?^ = -C2 = ^, (4.20) 



4.3 Four loops 

A novel feature that we first encounter at four loops is that the correlation function F^^^ receives 
non-planar corrections 

Fi*>=F^% + ±Xi,. (4,21) 

c 

The two functions -Fg=o and -Fg=i have the same general form (12. 5p and (12. 6p and are defined by 
the polynomial 

P^'\x„...,xs) = P^% + ^,P^^,. (4.22) 

c 

By construction, the polynomials -Pg=o and Pg=i are invariant under S's— permutations of the four 
external points xi, . . . , X4 and the four integration points X5, . . . ,xs- 

Similarly to three loops, the explicit form of the planar polynomial Pg=Q can be found from 
the requirement for the correlation function in the planar sector to have the correct asymptotic 
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behaviour in the hght-cone hmit 3^12; 2:23, 3^34, a^li — > 0. The result is [2] 

0(4) 12222222226 
-T g=0 — 24"^12-^13"''16-^23'^25-^34-^45'^46-^56-^78 

,12222222244 
-|- gJ-i2-^l3J'l6"''24-^27"''34"''38-^45-^56'^78 

1222222222222,c //i 9q\ 

~ 16'^12-^15'^18'^23'^26-^34'^37-^45'^48-^^ ~r '-'8 permULaXlOnS . (^4.Zoj 

In the non-planar sector, the same requirement turns out to be less restrictive but it allowed us 
to determine the non-planar polynomial Pg=i up to four arbitrary constants which are expected 
to take rational values. In order to fix these constants, one would need more detailed information 
about the properties of the correlation function. To save space, here we do not present the general 
expression for -Pg=i, it can be found in Ref. [2]. 

In fact, for our purposes we only need the expressions for the polynomials Pg=Q and Pg=i in 
the short-distance limit X2 Xi and 0:4 — )■ X3, Eq. fl4.4p . For the planar polynomial we find from 
( I4.23P for X2 = Xi and X4 = X3 

0(4) 8/222244, 222262i/i2222242 

-< g=0 ~ ■^13l,"''15-^18"''36-^37'^58"*'67 "''16"''17"''36"*'37"''58"''67 ^"''17"''18"''36"''37"''56"''58"''67 

i/i2222242io2222422|/|22222222 
"T ^'''16"''17'''37"''38"''56"''58"''67 ^ ^■^17-*'18-^37"''38"''56'^58"''67 ^"''17"''18-^36"''38-^56-^57-*'58-*'67 
io2 2 2 2 2 2 2 2 8\io6/'2 4 4 2 6 io4 2 2 4 4 2 

_l_ / nn ry> rf* /y <-y> <^ <-y> <^ <-y>^ I _l_ /If \ ly ly y ly ly _]_ J ly ly /y ly ly /y 

"T ^"''16"''18"''35"''37"''57"''58-^67"''68-^13j ^-^13 V"''17"''18-^37"''38"''56 ^"''17"''18-^36'^38-^56'^57 
io2442 42ir)2442 2 2 2io24242 2 2\ 

_l_ / ry ry ry /y /y ly —1— / If If If If If If If —X— / If If If If If If T* I 

-t- ZX;^gX;^gJ.37X3gX5gJ.57 -t- Z.Xi-jd,i^X^QX2,jX^QX^'jX^^ -f- Z,Xi^Xi^X2,QX2,j-L^Q-L^^d'QT ) 
io4/262244| 224624, r,2242424 

_l_ / ly I ly y 'T* If nf If —I— If nf If ry 'f —I— / If If If 'f 'f 

"T ^"''13i,"''17"''18'^35-^36'^37"''56 "''15"''16"''18"''37-^38"''56 ^'''16-^17"''18-^35-^37"''38"''56 
io2242622ir)2424422io2624222 

_l_ J /y ry ■T* nf ly if If —I— /if If If If If If If _l_ /If 'f If If If If If 

-f Z^^i^X^qXi-jX^qX^^X^qX^-j -t- jLXi^XirjXi^X2,QXy^^X^QX^'j -t- Z,XijXi^X2,^X2,qX2,tX^qX^'j 
io4424222 22442222\|^2/2262262 

_l_ / ly ly If If If nf If If nf If If If If 'y nf l —I— /I if I if if ly if if ly If 

"T ^•^17-^18"''35"''36-^38-^56-^57 "''16'^17-^18-^35-^36-*'37"''58"''67^ ^-^13 l,"''15"''16"''18-^35-^36'^37"''56 
io2264242io2442622io4244422 

_l_ / ly y /y " /y /y ly /y _l_ / ly /y ly /y ly ^ rf /y _[_ / /y /y <-y» /y /y rf /y 

-|- ^"f-i6-^17"''18"''35-^36"''37"''56 ^ ^"''15-^16-^18"''35"''37"''38'''56 ^"''16"''17"''18"''35"''37"''38"''56 

- xl^xl^xl^xl^xl^xl^^xl-jxl^xlf^) + ^4 permutation . (4.24) 

Here the 6*4 permutations are needed to restore the symmetry of -Pg=o under the exchange of the 

integration points X5, . . . ,X8. The polynomial -Pg=o is also invariant under the exchange of the 
external points Xi and X3, so that it has an S2 x 5*4 permutation symmetry. Notice that the 
relatively long expression for -Pg=0) ^ compared with (I4.23p . is an artifact of decomposing the 5*8 
permutations into S'4 x 5*4 ones. 

For the non-planar polynomial the situation is just the opposite. The expression for -Pg=i is 
very lengthy whereas in the short-distance limit it takes a remarkably simple form. Namely, the 
four different polynomials that accompany the four arbitrary constants in the expression for -Pg=i 
become proportional to each other at X2 = Xi and ^4 = ^3, leading to [2] 

-Pg=l — J/Sa "^§=1 ~ 6*"^^^ ^13(^56^78 + ^57^68 + ^8^67) Y\. ^l^'^Si ? (4.25) 
^3^^4 j=5,6,7,8 

with Cgt^]^ given by a linear combination of the four rational constants mentioned above. 

Then, we substitute the polynomial P*^^-* into (14. 3p and use (14. 2p to obtain the expression for 
Going to the corresponding integrand X4, Eq. (14. 6p . we find 

X4 = X4,g=o + -^^4,,=i . (4.26) 
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Here, the expression for the non-planar correction reads 



4! (4vr2)4 (a:f,a:|5) . . . {xi.xj,) U 



5<i<j<S ^ij 

^ ^ ^8 / ^2 ^2 1^ ^2 ^2 1^ ^2 ^2 '\ 
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A \ ( A 'rr2 ^ 4/^^2^2\ /'™,2^2\ ^2^2^2^2^2^2 ' ^ ^ 

^. l^^/l ; l-'-lS-^SSJ • • • l-^lS-^SSJ "''56"''57"''58-^67-^68"''78 

where in the second relation we replaced -Pg=i by its explicit expression fl4.25p . In the planar 
sector, the analogous expression for Xi,g=o is much longer since it involves the P— polynomials at 
lower loops: 



6 xf 

-4,g=0 



-^5,6,7,8 4 -^5,6,7 



/I I /y»2 ^2 /j-»2 /vi2 fy»2 /j>2 13 ^2 ^2 fy»2 

^. J'56"''57't'58-^67"''68'''78 -^bG-^bT^Gr 



+ 5*4 permutations (4.28) 



_3 4 P5,6-P7,8 , 1o^™4^2-P5,6 p:^^_4n3 
4 

where P5,6,7,8 = Pg=o given by fl4.24p and the right-hand side is symmetrized with respect to 
all 5*4 permutations of the integration points X5, . . . ,xs- Replacing the P— polynomials by their 
explicit expressions, we obtain a lengthy result for X4^g=o- 

At the next step, we restrict all the integration points X5, . . . , xs to the vicinity of the external 
point Xi and simplify the integrand X4 by replacing xf^ — )■ xf^ (with i = 5, . . . , 8). Denoting the 
resulting function X4, we find from f l4.27p and fl4.28p 

/^^ ^2 I ,7i2 ^2 I ^2 ^2 

2; _ g=l •^56-^78 "T •^57-^68 ~r ■^58-^67 29) 

yi I / /l'r]-2^4 /v>2 ^2 /v>2 ^2 /vi2 ^2 /vi2 ^2 /ji2 fy»2 ^ ^ 

'i. l^'i/i ; •^l5-^l6-^17'^18-^56-^57-^58-^67-^68'^78 

in the non-planar sector, and 

^ 

A] { /1_'7t2 ^ 4 /-v*^ <T^4 ^4 /y»4 /j>2 fy»2 /ji2 ^2 /vi2 fy»2 
^•l^" J •^15-^16-^17'''18'''56"''57'''58"''67"''68-^78 
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[^•^17-^ 56-^ 18 "T ^"''15"''16"''56"''18 °-^16"''17"''56"''18 ^"''17"''56"''57-^18 ^"*'17"''56"''l 



A 

as 
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10 2224 i^4422|^24222|^22422 

IZXggXgyXggXgyXlg "I" ^X j^yXggXg^X ]^g "I" ^X j^gX XggXgyX ]^g "I" "±X j^yXggXggXgy X ]^g 

ir>24222 on 222222 Ar 2 2 2 2 2 2 
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17-^56-^58-^67-^18 ^"-'■^17-^56-^57-^58-^67-^18 ^'^■^16-^57-^58-^67'^68-^18 
222222i/i22422i 2262i/i22242 

/y /y /y /y /y _l_ /\_'y /y /y /y rf _l_ /y /y ly " /y _l_ /I /"vt /y> " /y'-' rf rf 

56-^57"''58-^67"*'68"''18 ^"''15"''16-^17-^56'^57 "''16"'' 17-^ 58-^67 ^"''16"''17"''56"''58"''67 



— 12x^gX^7X5gX57X5gXg7 — 21x5gX5yX5gXgyXggX7g] + 6*4 permutation (4.30) 

in the planar sector. We verify that the functions X4 g=o and X4 g=i transform covariantly under 
conformal transformations and do not depend on the external point X3. 

Then, to find the four-loop correction to the Konishi anomalous dimension fl3.33p . or equiva- 
lently the coefficient C3, we substitute X4 into fl3.30p and obtain: 

C, = ^^1^ J d'-'^xed'-'^x,d'-'^xs (x4,g=o + ]^^4.g=i) . (4.31) 
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By construction, C3 is dimensionless and is expected to be finite as e — )■ 0. We do not remove 
the regularization, however, since it is more advantageous to expand the right-hand side of (14.311) 
into a sum of basis three-loop (divergent) integrals and evaluate each of them separately. 

Let us start with the non-planar contribution to (14.311) . Setting for simplicity xl^ = 1, we 
find from fOOD . (|43T|) 



J4) 



a 



3,g=l 



10247r6 



j4-2e 



Xs 



d^-''xed^-''x7d 



2222222 



(4.32) 



The integral on the right-hand side of this relation is shown diagrammatically in Fig. El It is 
finite at e = and, most importantly, it corresponds to a planar graph. This fact allows us to 
introduce dual momenta and represent the integral in the form of a planar dual (momentum) 
graph, as shown in Fig. |3l The main advantage of using the dual representation is that the 




X8 




Figure 3: Diagrammatic representation of the integral (I4.32p in the representation (left) and 
in the dual momentum representation (right). In what follows, all momentum integrals are shown 
in blue. 

resulting three-loop momentum integral can be easily evaluated using the MINCER package [3T] . 
leading to 



C, 



1024 



X (2OC5 + 0(e)) . 



(4.33) 



We recall that the constant c"^li is expected to take rational values. 

Let us now examine the three-loop integrals in fl4.3ip generated by the integrand X4g=o, 
Eq. fl4.3Up . We shall denote the corresponding contribution to the right-hand side of (14.311) by 
0*3,5=0- The three- loop propagator integrals appearing in the calculation belong to the following 
family of three- fold integrals, with various integer indices Oi, . . . , Og 



0(0.1, . . . , ag) 



d^-''x^d^-''x7d 



.2 ^ 



.^^r^{xiY^{x 



[xlY^'ixl^fHxt, 



)as(^2^)a, ' 



(4.34) 



where we put X5 = for simplicity. Notice that the indices can take both positive and negative 
values. In the latter case, the corresponding term appears in the numerator. For arbitrary choices 
of the indices some of the integrals (I4.34p may be non-planar and so can not be rewritten as 
dual momentum integrals of propagator type. However, most importantly, in our case all the 
integrals appearing in (I4.3ip are planar and thus of propagator type. 



21 




^ -Q- <xx> 




Li Pi P2 Pz P4 

Figure 4: The dual momentum master integrals defining the four-loop Konishi anomalous di- 
mension. 

To evaluate each of them we apply the standard technique called integration by parts (IBP) 
[33] (see Chapter 5 of |3l] for a review of the method) which provides the possibility of representing 
a given integral of this family as a linear combination of so-called master integrals. We found 
that the calculation of C3^g=o involves only 5 master integrals, all corresponding to planar graphs. 
Therefore, introducing the dual momenta ki = Xi — Xj+i, we can rewrite these integrals as the 
dual momentum-space integrals Li, Pi, . . . , P4 shown in Fig. |H 

The resulting expression for C3,g=o is 

/ 9 3 45 165e 16A7e^ ^ ^ . 4. 

C3 g-o = Li \ \ 54e^ + O (e^) 
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The three-loop massless propagator master integrals in this relation were evaluated many years 
ago [33]. Using the known results [27] for the master three-loop integrals Li, Pi, . . . , P4 we find 
from (05]) 

39 9 45 

C3,g=o = y- 4(3 + ^(5 + 0(e). (4.36) 

Substituting f l4.33p and f l4.36p into fl3.33p we finally obtain the following result for the four-loop 
correction to the Konishi anomalous dimension 

(4) 39 9^ 45^ r ^ 
7g^ = -^ + 4C3-yC5 + ^C5, (4.37) 

where r = — 5Cg'2i/256 is an undetermined rational constant. The planar correction to 7^"* is 
in agreement with the known result [3 O 13 |15l [161 HI]- The non-planar correction has been 



^°A11 the integrals appearing in (|4.3ip can be handled by MINCER [3T]. In addition to the single integral ()4.32p 
contributing to the non-planar correction to the anomalous dimension, we found 76 integrals in the planar sector. 
Evaluating them with the help of MINCER, we arrived at the same result ()4.36p . However, we would like to 
emphasize that at five loops the reduction of more than 17000 integrals to master integrals can only be done by 
a direct application of the IBP method (see Sect. [5]). 
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computed in [26], and the result confirms our prediction. Moreover, it allows us to fix the value 
of the unknown constant in fl4.37p : 

r = -1^ . (4.38) 

In this section we have demonstrated the high efficiency of our method for computing the 
Konishi anomalous dimension at four loops. We emphasize once again that we do not use the 
conventional Feynman diagram technique. For comparison, the direct calculation of [13 ITB] 
involves hundreds of = 1 super-graphs, each giving rise to a number of Feynman integrals; in 
the calculation of [T7] the number of contributing Feynman graphs exceeds 130 000. 



5 Konishi anomalous dimension at five loops 



It is straightforward to extend our analysis to five loops. In this case, the correlation function 
F^^^ has the following form 



p(5)( 



<5) 



1 



I J p(^) J P 



1 



(5) 



(5.1) 



In what follows, we shall restrict the discussion to the planar sector only. As was shown in 
Ref. [2], the five- loop correlation function in the planar limit is given by 

d'^X^d^XQd'^Xjd^Xsd'^Xg P^^\xi, . . . , Xg) 
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(5.2) 



5! (-47r2)5 

where the polynomial P^^^ is invariant under the 5*9 permutations of the four external points 
xi, . . . , 0:4 and the five integration points x^, . . . ,xg. It is given by the following expression: 

d(5) 122224222222222222 
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j^2^13^16^17^19^26'^27-^28-^29-^35-^38-^39-^45'^46-^47'^49-^57-^58-^68 ^ ^9 pcrmutatioUS , (5.3) 

where the relative coefficients follow from the requirement for the correlation function to have 
the correct asymptotic behaviour at short distances. 

The analysis goes along the same lines as at four loops. We start with examining the integrand 
(15.21) in the double short-distance limit X2 — )■ Xi and X4 — )■ X3 and, then, apply f l4.2p to identify 
the five- loop integrand of I^^^: 
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(5.4) 
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where ^5,6,7,8,9 is the polynomial P^^^ evaluated at X2 = xi and X4 = x^ and the remaining 
P— polynomials were defined earlier (see Eqs. (14.71) . (14.121) and (I4.28P ). The expression on the 
right-hand side of (15. 4p is symmetrized with respect to S5 permutations of the integration points 
X5, . . . ,xg. To save space, here we do not present the explicit expression for X5. 

At the next step, we simplify the expression for X5 by choosing all integration points x^, . . . ,xg 
to lie in the vicinity of the point xi. This is equivalent to sending the external point to infinity 

— )■ cxD with all remaining points fixed 



%= lim X5(xi,a;3; X5, . . . ,0:9) . (5.5) 
0:3— ^00 

Finally, we apply (I3.33P to express the coefficient C4 as the following four-loop integral 

2tt^ f _ _ _ _ - 
= . 2 N9^A. / "^'xed"^ '^'xj(f ^'xsd'^ ^'xgX^ixi; X5, . . . , Xg) . (5.6) 



As before, to simplify the calculation we put X5 = 0. 

Replacing I5 in (15. 6 p by its explicit expression, we find that C4 is given by the sum of more 
than 17000 four-loop two-point Feynman integrals. All of them belong to the following family of 
four-fold integrals, with various integer (positive and negative) indices ai, . . . , 014 
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As in the four-loop case, to evaluate each of them we apply the IBP method [33]. To solve the IBP 
relations, i.e. to represent every integral on the right-hand side of (15. 6 p as a linear combination 
of master integrals, we apply the C++ version of the code FIRE [35]. In this way, we found that 
C4 is given by a linear combination of 22 master integrals: 



= W44M44 + WgiMgi + WsqMsq + W31M31 + Ws^M^ir^ + W22M22 + W32M32 

+ U733M33 + W34,Mu + W25M25 + W23M23 + W27M27 + W24M2A + ^26^26 

+ WqiMqi + W21M21 + W12M12 + WiiMii + W14M14 + W13M13 + Wih + W2I2 , 



(5.^ 



with the coefficient functions Wi defined below in Eq. (15. 9p . Among the master integrals only 
two, Ji and I21 are associated with non-planar graphs (see Eqs. (IB.ip below). The remaining 
20 master integrals M44, . . . , M13 correspond to planar graphs. This allows us to introduce the 
dual momenta ki = Xi — Xj+i and represent the same integrals as four-loop propagator master 
(momentum) integrals shown in Fig. |5l The latter i nteg rals were calculated recently in [27] as 
an e expansion up to transcendentality weight seven, [jj The explicit expressions for the planar 
integrals M44, . . . , M13 can be found in [21]. To save space, we do not present them here. 
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At the moment, results for the master integrals are known up to transcendentality weight twelve [36j . 
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The corresponding coefficient functions are given by 
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Here the series expansion of Wi is truncated at the order in e which is related to the maximal power 
of 1/e in the expression for the corresponding basis integral Mj in fl5.8p . so that the right-hand 
side of (15 -Sp can be evaluated at order 0(e°). 

The two non-planar master integrals Ji and I2 entering the right-hand side of (15. 8p are shown 
diagrammatically in Fig. [H] and their explicit form can be found in (IB.ip . These integrals are 
evaluated in Appendix [B] leading to 

, 5 _g 1 o 2 , or- A , / 13 4^ r^i ^ 2 , ir^r^ ^ _2 ^ , ^45 
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Their coefficient functions in (15. 8p are 
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Figure 5: Diagrammatic representation of the planar basis integrals in the dual momentum 
representation. Blue line denote momentum propagators with the momentum k = Xi — xj. 



Finally, we combine together the relations ( 15. 9p - (15. lip , make use of the results of Ref. 
for the master integrals shown in Fig. [5] and obtain from (15. 8p the following result for C4 or 
equivalently, five-loop Konishi anomalous dimension 



(5) 



237 27 



135 



945 



16 



(5.12) 



This relation is the main result of the paper. It is in perfect agreement with the prediction of 
the integrable models |5l El [71 E] ■ 

6 Conclusions 



In this paper we have developed a new efficient method for the computation of the Konishi anoma- 
lous dimension at higher loops. It does not use the conventional Feynman diagram technique 
with the associated very large number of contributing graphs and Feynman integrals. Instead, we 
exploited the recently discovered new symmetry of the four-point correlation function of A/" = 4 
SYM stress-tensor multiplets to predict the form of its integrand as a linear combination of a 
small number of relevant diagrams. Then, we examined the asymptotic behaviour of the loga- 
rithm of the four-point correlation function in the double short- distance limit and related the 
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Figure 6: Diagrammatic representation of the non-planar basis integrals Ji, I2 and two auxiliary 
integrals /3(0), /4(0) defined in Eqs. fIB.ip and fIB.lOp . respectively. The fine with the index 2 
denotes a square of scalar propagator while all the remaining lines stand for l/x"^. The 

points Xl and are external and integration goes over the points Xq, . . . ,Xg. 

Konishi anomalous dimension to its leading logarithmic singularity. Finally, by analyzing the 
expected singularity of the logarithm of the correlation function in this limit, we were able to 
lower the loop order of the contributing Feynman integrals by one, that is to express the Konishi 
anomalous dimension at i loops in terms of finite two-point integrals at {i — 1) loops. Going 
through these steps, we obtained the five-loop Konishi anomalous dimension in the planar limit 
as a sum of 22 master four-loop two-point integrals. Replacing the master integrals by their 
explicit expressions we arrived at an analytic result for this anomalous dimension which agrees 
with the integrability prediction |S1 El EJ [S] . 

At present, the expression for the integrand of the four-point correlation function is known up 
to six loops in the planar limit [2]. In our calculation of the Konishi anomalous dimension we made 
use of this expression to five loops only. By applying the method developed in this paper, it is 
straightforward to extend the analysis to six loops and to express the six-loop Konishi anomalous 
dimension in terms of five-loop two-point integrals. The evaluation of such integrals is still an 
open problem, not because of their number, but because the IBP reduction to master integrals is 
a very complicated problem at this level. Still, we are optimistic that further development in this 
direction will eventually make the six-loop calculation possible. Likewise, no prediction for the 
six-loop Konishi anomalous dimension is available from AdS/CFT considerations, and to obtain 
it using the existing integrability approaches appears to be a rather non-trivial task. 

Another result of our study is the prediction of the non-planar correction to the Konishi 
anomalous dimension at four loops in the form r^^/N^ with r being an undetermined rational 
number. Our prediction is in full agreement with the result of the direct Feynman diagram 
calculation in [26], which also allows us to fix the value of r = —135/2. It is interesting to note 
that in our approach the non-planar correction at four loops originates from just a single and 
very simple three-loop propagator integral shown in Fig. [3J We would like to emphasize that the 
non-planar 0{1/N^) correction to the four-point correlation function derived in Ref. [2] depends 
on four arbitrary rational constants. The parameter r is given by a particular linear combination 
of these coefficients. To fix each of the four coefficients we need three more relations. They 
can be obtained from the comparison of the non-planar corrections to the twist-two anomalous 
dimensions computed in Ref. [M] with the analogous results from the OPE analysis of the non- 
planar correction to the four-point correlation function. As explained in Ref. [2], the perturbative 
corrections to the correlation function have an iterative structure at higher loops. In application 
to the Konishi operator, this implies that the non-planar 0(1/ N^) correction to its anomalous 
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dimension is uniquely defined at all loops by the values of these four coefficients. We would like 
to mention that starting from five loops, the anomalous dimension receives 0{1/N^) corrections. 
The method proposed in this paper can be equally applied to the study of such corrections. 

The Konishi operator is just the first in an infinite series of twist-two operators, all appearing 
in the OPE of two Af = 4 SYM stress-tensor multiplets. The four-point function that we use for 
the evaluation of 7ac contains the information about the whole spectrum of anomalous dimensions 
of twist-two (as well as higher twist) operators. However, in order to extract it from the OPE, one 
needs to either evaluate analytically all relevant higher-loop four-point conformal integrals, or at 
least to work out their asymptotic expansion in the double short-distance limit. This problem is 
not yet solved in full generality beyond two loopS but it undoubtedly deserves further attention. 
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Appendix A IR rearrangement in coordinate space 

In this appendix we explain in detail the method that we employ in our calculation of the Konishi 
anomalous dimension. It represents an extension of the so-called infrared rearrangement method 
(IRR) [38j to the coordinate space. 

To describe the method, let consider as an example the following four-loop integral in Eu- 
clidean D— dimensional space-time (with D = 4 — 2e) 

^K-^IZ) 2D I -2 ,,2 .,2 ,,2 .,2 ,,2 .,2 ,,2 ,,2 ,,2 „,2 ,,2 ' l^'-^J 



.-TT-^j_y f rft^ <-y»-" <-y»-" <-y»-" /y'-' rr>^ ry>^ /y*" nf'-' /y- 

n J Xi5X;^gXi7X;^gX35X3gX37X3gX5gXggX78X57 



We would like to stress that Xi are true coordinates in the configuration space and, therefore, 
/(xis) is different from the conventional integrals that one encounters in the dimensional regu- 
larization in which case all distances in the denominator appear in the power (1 — e). 
The integral (]A.1|) has a simple pole in e 



/(Xis) = (x 



-4e 



13; 



- + 0(e°) 
e 



(A.2) 



It comes from the integration over the region where xs, . . . , xs are all close to xi and from the 
symmetrical region where 0:5, . . . ,X8 are all close to X3. Since the integration variables as true 



^^Although at three loops, the inverse process has been partiaUy done, namely a prediction has been made 
for the three-loop correlation function in the limit u but with finite v 37J, by making use of the three- loop 
twist-two, arbitrary spin anomalous dimensions predicted in |24) together with known lower-loop twist-two data, 
and a conformal partial wave analysis. See appendix [C] for more on higher spin twist-two anomalous dimensions 
and the four-point correlation function. 
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coordinates in the Euclidean space, the pole 1/e has to be qualified as an UV divergence. Notice 
that the integrand of /(xia) coincides (up to an overall normalization factor) with f l4.27p . As 
a result, the residue C defines the four- loop non-planar correction to the Konishi anomalous 
dimension. 

In general, the UV divergences in the coordinate space come from regions where the integrand 
considered as a generalized function of Xi (tempered distribution, i.e. linear functional on a space 
of test functions) is ill-defined. In our example, the product of x^— factors in the denominator 
of (lA.ip turns out to be unintegrable in a vicinity of the two external points, Xi and x^,. In the 
first case, we consider the product 

F{xi, X5, . . . , Xg) p K p K p K K (A. 3) 

fT»-" <-y»-" rtf-^ ly" ry" ly" ry^ ly^ 

1 5 1 6 1 7 1 8 5 6 68 78 5 7 

as a tempered distribution. Its divergent part is described by an UV counter-term E2l 

C 

A(xi,X5, . . . ,X8) = i^^K^i - ^5) ■ --Si^xi - xg) , (A.4) 

with the constant C determined below. Similar counter-term A(x3, X5, . . . , xg) describes singular 
behaviour of the integrand (lA.ip in the vicinity of X3. Thus, the pole part of (lA.ip is just twice 



the factor C/(2e) in f EIij) 

1 = j rf^xs . . . rf^xg [A(xi, X5, . . . , xg) + A(x3, X5, . . . , xg)] + 0(e°) = ^ + 0(e°) , (A.5) 



leading to flA.2p 



To evaluate the constant C in (IA.4P we apply the infrared rearrangement (IRR) method orig- 



inally proposed by Vladimirov in Ref. [38] in the momentum space. It makes use of the fact that, 
for an infrared finite but logarithmically UV-divergent Feynman integral without sub divergences, 
the contribution of the counter-term is just a constant. The idea of IRR is to set the external 
momenta to zero and then, in order to avoid the appearance of IR divergences, to introduce an 
external momentum (or a mass) in such a way that the calculation becomes simpler!^ 

Applying the IRR method to ( lA.ip . we should have transformed the integral /(xia) to the 



momentum space via Fourier transform. However we will not do this for the following two 
reasons. First, the resulting momentum integral will be four-loop one while we can obtain 
the residue C from three-fold integral only as described below. It is well known |20] that the 
evaluation of the UV pole part of a given £— loop momentum-space Feynman integral can be 
reduced to evaluating massless propagator (£— 1)— loop Feynman integrals to order e°. However, 
as was already mentioned, the integral ( lA.ip is different from the conventional Feynman integral. 
In particular, the 1/x^ factors on the right-hand side of ( lA.ip are replaced in the momentum 



representation by the factors of l/(fc^)^~'^ depending on e. These are much more complicated 
objects, both from the point of view of an IBP reduction and evaluating master integrals, so that 
it is the second reason why we want to stay in coordinate space. 

Let us apply the IRR method to ( IA.3P in the coordinate space and treat the coordinates Xi, X5 
as external and X6,X7,Xg as internal points. Notice that setting an external momentum to zero 



^■^In this example, A contains a simple pole in e. In a general situation, this would be a finite linear combination 
of negative powers of e. 

^^If it is not possible to avoid such IR divergences one can remove them immediately by the so-called R*- 
operation |391 120] but we do not meet such a complication in our calculations. 



29 



corresponds to integrating over the corresponding coordinate. Then, the constant C in flA.4|) can 
be obtained by integrating the both sides of flA.Sp with respect to internal points 



F{xi,X5) 



jD 



22222222 

/->'■ ryi^ rf'-^ rr*^ /y^ ry^ rr*^ 



^5(xi-X5) + 0(e°). 
2e 



(A.6) 



The integral on the left-hand side depends on the two external points and is of propagator type. 
We can check it has no IR divergences, i.e. divergences at large values of coordinates, and has 
the following form by dimensional arguments 



F(xi,a;5) = /(e)- 



1 



(A.7) 



Here the only source of the simple pole in e is hidden in the second factor (which is considered 
as a distribution) so that /(e) is analytic in a vicinity of the point e = 0. The simplest way to 
reveal the 1/e pole of the distribution l/{x\^y'^^'^ is to take its D-dimensional Fourier transform 
with a help of the identity 



1 




_(X2)\ 





d^x e 



ipx 



_ 4^/2-^r(D/2- A) 

(X2)~ ~ T(A) (p2)13/2-A • 



(A.^ 



In particular, for A = 2 + 3e we find from (I A. 71) (for xs = 0) 

4-4T(-4e) 1 /(O) 



^[F(xi,0)] = /(e) 



r(2 + 3e) (p2)-4e 



4e 



0(e° 



(A.9) 



At the same time, replacing -F(xi, 0) by its expression (1A.6P we obtain the left-hand side of this 
relation as C/(2e) + 0(e°) leading to 



C 



:/(0) 



-F(xi,0) 



(A.IO) 



=1,Z)=4 



It is easy to see that the integral F{xi,x^), Eq. (IA.6|) . corresponds to a planar graph shown in 
Fig. 131 After going to the dual momenta, we find that it coincides with the graph N2 of Baikov 
and Chetyrkin [27] leading to 



C = -10C(5). 



(A.ll) 



Appendix B Non-planar master integrals 

In this appendix, we evaluate the two non-planar master Euclidean integrals f IS.lOp . They have 
the following form (with D = A — 2e and x\ = 1) 

J _ e^ f d^xed^xjd^xsd^xg _ Qi ^ ^ ^ce + 0(e^) 

/I J XiqXiqXq-jXq^XjX-j()X^X^^ t 



^2D 



d^Xed^Xjd^Xgd^Xg '^2 , , , , 2 



'T'2 { rp'^ '\2,-yj2 ,-yi2 ^2^2 ^2^2 ^ 



+ b2 + C2e + 0{e'). (B.l) 
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Here we introduced the factor in front of the integrals to avoid the appearance of terms propor- 
tional to In TT and Euler's constant 7 in the right-hand side. The diagrammatic representation of 
Ii and I2 is shown in Fig. Ei Both integrals develop poles 1/e but their origin is different. For 
the integral Ji it comes from integration over xq,Xt,Xs,Xq going to infinity simultaneously and, 
therefore, has an IR origin. For the integral I2 the pole comes from integration at short distances 
Xig — )• and has UV origin. 

Substituting Eqs. (15. 9 p - (15. lip and ( IB.ip into (15. 8 p and making use of the results of Ref . [27] 
we finally obtain the following expression for C4 
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(B.2) 



Here the constants a^, hi and Cj describe the contribution of the two non-planar master integrals, 
Eq. fIB.ip . We recall that C4 defines the five- loop correction to the Konishi anomalous dimension 
and, therefore, it should be finite for e — )■ 0. The condition for the 1/e^ and 1/e poles to cancel 
inside C4 leads to two relations between the coefficients and 6j. As we shall see in a moment, 
these relations are indeed satisfied. 

Let us start with the leading 0(l/e) term on the right-hand side of fIB.ip . The simplest way 
to compute the residue at the pole is to Fourier transform the integral into momentum with the 
help of f lXSjl . Notice that the expressions on the right-hand side of (IB.ip are valid for x\ = 1, 
but their dependence on x\ can easily be restored from dimension analysis. In this way, we find 
from the first relation in (IB.ip 



^[/i]=^ ^(x?)-^^ + 0(6°) ={Ua, + 0{e)){p 



L e 



-2+5e 



(B.3) 



This relation implies that the coefficient ai can be obtained from the Fourier transformed integral 
J-" [/i] evaluated at D = 4 dimensions. Transforming the integral Ji into the momentum represen- 
tation we find that J-" [Ji] coincides (up to a factor of 16) with the conventional four- dimensional 
momentum Feynman integral denoted A^^o in Ref. [27] 



-^[/i] = 16 



Comparing this relation with (1B.3P we find 




16 X (20C5 + O(e))(p- 



2\-2+5e 



(B.4) 



ai = 5(5 • 



(B.5) 



^^Since the Euclidean integrals in (IB.1[) are positive definite, this explain why their residues at the pole, ai and 
02, have opposite signs (see Eqs. (jB.5l) and (IB.9|) below). 
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Let us now turn to the integral I2 in flB.l|) and Fourier transform it 



e 



4(f + 0(e; 



-l+5e 



(B.6) 



To identify the momentum integral corresponding to J^[l2] we have to Fourier transform all 
factors in the denominator of I2 including l/(x^g)^. In that case, we find from ( ]A.8|) 
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The fact that the residue at the pole in this relation does not depend on the momentum p implies 
that the corresponding line in the Feynman diagram shrinks to a point. As a result, 
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e 





-x(20C5 + O(e))(p^ 
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Here in the second relation we redrew the same diagram, so that it takes the form of the diagram 
N2 in the notation of Ref. [27j. Comparing the last relation with flB.6p we conclude that 



a2 = -20C5. (B.9) 

Let us now compute the subleading terms in the expansion (IB.ip . To this end, we introduce 
the following auxiliary integrals (with D = A — 2e) 
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16 17-^ 19-^ 67-^ 68-^ 79-^ 7-^8-^ 89 

with K, being a parameter. The diagrammatic representation for these integrals (for k = 0) is 
shown in Fig. [6l For finite k the two integrals are finite as e — )■ and, therefore, they admit an 
expansion in powers of e and ne 



lAn) = bi + e(ci + ndi) + 0(e^ 



3,4) 



(B.ll) 



where we put xl = 1. Notice that the leading 0(e°) term does not depend on k, whereas the 
0(e) term is a linear of function of n. 

A distinguishing feature of the integrals ( IB.lOp is that, in the special case k = 0, the IBP 
relations allow us express /3(0) and /4(0) in terms of two master integrals Ii and I2, Eqs. (IB.ip . 
More precisely. 
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Once the auxiliary integrals (IB.11|) have been computed, we could use these relations to obtain 



the needed coefficients bi, 62 and Ci, C2. 

One may wonder why we introduced the parameter k into the definition of the integrals fIB.lOp 
if we only need its value at k = 0. The reason for this is that, as we will see in a moment, it is 
much easier to compute the integrals (1B.10|) for the two special values k = 1/2 and k, = 1. Then, 



taking into account that /^(k) is a linear function of k at order 0(e), Eq. (]B.11|) . we find 

hiO) = 2h{l) - hil/2) + 0{^) = bi + ec, + 0(e2) . (B.13) 

In what follows, we shall evaluate /«(!) and 7,(1/2) and, then, apply this relation to compute 
63, 64 and C3, C4. 

Let us consider the integrals fIB.lOp for k = 1. In this case it is easy to see that the integrand 
is given by a product of factors l/{x'^Y~'^ which coincide with scalar propagators in D = 4 — 2e 
dimensions. As a result, upon the Fourier transform, the integrals ^-"[73(1)] and ^-"[74(1)] are given 
by conventional four-loop momentum Feynman integrals. In this way, we find that the integral 
^-"[73(1)] coincides with the master integral M45 in the notation of [27] 



73(1) = G^M45 = 36C3' + 6(108(3(4 + 288C3' - 378 Ct) + 0{e') , (B.14) 

where the additional factor Go = e^^ r(l + e)r(l - e)Vr(2 - 2e) = 1 + 2e + 0(e2) is inserted to 
convert the result of Ref. [27] obtained in the G— scheme to the regularization scheme used in 
(IB.lOp . The second momentum integral ^-"[74(1)] is not a master integral. We applied FIRE to 
reduce it to the master integrals of Ref. [27] and arrived at the following result: 



(3 - 4e)(l - 4e)(4 - 5e)(3 - 5e) (10 - 105e + 326e2 - 319e3) 
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4(1 - 26)(2 - 36)(1 - 36)(3 - 46)(1 - Ae) 2(1 - 46)(1 - Se) 

3(l-e)e4 3(1 -e)e 

_ (1 - 2e)(2 - 36)(3 - 56)(2 - 96)(7 - 19e) 1 - 5e 



+ M.^^-^-^(^-^:^?-^^^^ + mJ-^^^^^^ . (B.15) 
[l-e)e^ 3(l-e)e3 ^ ' 

Replacing the basis integrals by their explicit expressions we get 



74(1) = 36C3' + e (1O8C3C4 + IO8C3' + ^Ct) + Oie') , (B.16) 
where we put = 1. 

Let us now examine the integrals (IB.lOP for k, = 1/2. In this case, the special feature of 



the integral 73(1/2) is that the conformal weight of the integrand at the integration points Xj 
and xg equals the space-time dimension 4(1 — Ke) = D. As a consequence, performing inversion 
x^ —J- x^ /xf we obtain the following representation for 73(1/2) (at x\ = 1) 



J /-, ir)\ Q ^ C d Xgd X'jd XgC? Xg xi I XT/ \xs \ /-T) 

-'SU/ZJ — — ^ / (^.2 „2 „2 „2 „2 „2 „2„2„2 \l-e/2 ~ I \ 7 1 ' [^■^') 
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where on the right-hand side we depicted the corresponding Feynman diagram. Compared with 
the first relation in f lB.lOp . the product x^xf gets replaced here with x^x^. We observe that the 
diagram on the right-hand side of flB.17p contains a two-loop subgraph. As a consequence, the 
integration over Xj and Xs can be easily performed with the help of identity 



TT 



D 



d^xrd^^ 



(/y»2 2 /y»2 /j-»2 ^2 X 1 — c / 2 
'''67"''68'''78-^79"''89/ 



2 )l-e/'. 
69) 



[X, 



[6C3 + (9C4 + 12C3)e + 0(e^ 



(B.18) 



Substituting this relation into (IB.17P we find that the remaining integral over xq and Xg takes 
the same form as (IB. 180 leading to 



2M2 



h{l/2)= [6C3 + (9C4 + 12C3)e + 0(e2)] 



(B.19) 



It remains to determine /4(l/2). The corresponding Feynman diagram has the same form as the 
one for J4(l/2) (see Fig. [6]) with the only difference that each solid line carries the index (1 — e/2). 
It is easy to see that the integrals /4(l/2) and /3(l/2) look differently. Quite remarkably, as we 
will show later in this appendix, they coincide leading to 



J4(l/2) = /3(l/2) = 36C3' + e (IO8C3C4 + I44C3') + 0{e') . 

We would like to stress that this relation is exact and it holds for arbitrary e. 

Then, we combine the relations (IB.14p . (lB.16p . (IB.20p together and find from (IB. lip 

h{f^) = 36C3' + 6(108(3(4 + 288 /€C3' + (1 - 2k)378C7) + 0{e') , 
h{K) = 36C3' + e(l08 C3C4 + (180 - 72k)Cs^ - if (1 - 2k) Cr) + 0{e^) . 



(B.20) 



(B.21) 



Matching these expressions into (IB. lip we obtain the following relations for the coefficients 

63 = 64 = 86(3', 

C3 = 108(3(4 + 378(7, 



04 = 108(3(4 + 180(3^-^(7. 



(B.22) 



Their substitution into 
ci,C2 whose solution is 

ai = 5 (5 , 



+2|) yields a system of linear relations for the coefficients 61, &2 and 
5 



02 



Cl 



-20(5 
13 



b2 



378 
_ 10 
^189 



7r°- 13(3^ + 35(5, 



716-8(3^-40(5 



30 

4 



vr\3 - 91(3' + 195(5 - 77^ Cs + ^(7 + ^vr^ 



345 
T 



Co = VT (3 

15 ^ 



20 



16(3' - 80(5 + y vr2(5 - 520(7 



5 

54' 
20 



189 



-TC 



(B.23) 



Substituting these relations in ( IB.2P we verify the cancellation of poles in e and reproduce (I5.12p . 
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We complete this appendix with a proof of the relation /4(l/2) = 13(1/2). It relies on 
applying the cut-and-glue method of Ref. [331 I27j^ . Let us examine the integrals fIB.lOp for 
K = 1/2 + A/(10e) with A arbitrary. Dimensional analysis shows that the integrals have the 
following form 



1,(1/2 + A/(10e)) 



Ci(e, A) 



(3;2)l-./2-9A/10' 



3,4) 



(B.24) 



with Cj being some function of e and A. Then, the relation 14(1/2) = 13(1/2) implies that for 
arbitrary e 



C3(e,0) = C4(e,0). 



To show this, we consider the following Fourier integral 



J" 



1,(1/2 + A/(106)) 

(^2;2^1-<:/2-A/10 



c,(e. A) 

(Xf)2-^-' 



Q(e,A) ^ ''^^^\A pY' 



(B.25) 



(B.26) 



where in the second relation we applied ( IA.8I) . A crucial observation is that for small A the 
expression on the right-hand side develops a pole 1/A with the residue independent on the 
momentum p 



J" 



1,(1/2 + A/(10e)) 

(3,2)l-e/2-A/10 



A^ 



C,; e 



0) 



r(2 



+ 0(A°) 



(B.27) 



The integrals /,(l/2 + A/ (lOe)) (with i = 3,4) are described by the Feynman diagrams shown in 
Fig.[7]in the left column. All solid lines in these diagrams correspond to factors of l/(x^)^~'^/^~^/^°. 
Then, dividing J, by (xlY^^^'^^'^^^^ amounts to adding one additional line to the diagram con- 
necting the external points Xi and 0. The resulting diagrams are shown in the right column of 
Fig. [71 Notice that, up to changing the labels of the dots, these two graphs coincide (to see this, 
it suffices to rotate the lower diagram clock- wise by 271/3). After taking the Fourier transform, 
in the momentum representation, the momenta p and —p are injected into the external points 
Xi and 0, respectively, while for the remaining integrated points the corresponding momentum 
equals zero. As was argued in [27], the very fact that the leading asymptotic behaviour of the 
Fourier integral for A — ?■ does not depend on the momentum p implies that the residue at the 
pole 1/A is not sensitive to the choice of the external points. This is exactly what happens for 
the two graphs in the right column of Fig. [71 the only difference between these graphs is in the 
assignment of the external points. Since their residue at the pole 1/A are defined by the functions 
cs{e, 0) and C4(e, 0), we conclude that they are equal to each other leading to ( IB. 251) . 

To check numerically our analytic results for these two non-planar integrals we used the code 
FIESTA [41J which gave the precision of six digits. 



Appendix C Twist-two anomalous dimensions 

In Section !?^ we have applied the OPE (12. 8p to identify the contribution of the Konishi operator 
to the four-point correlation function in the double-short distance X2 — )■ xi and 0:4 — )■ X3. In this 

^^It turns out that this generahzed gluing is very close in its spirit to the strategy of ref. [ 40j where Feynman 
integrals were considered as distributions with respect to the parameter of analytic regularization. 
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Figure 7: Glue procedure for the integrals J3 and I4 defined in (lB.24p . All lines correspond to 
factors of with the same index a = 1 — e/2 — A/10. The dots with the labels Xi and 

describe the external points, the remaining four dots describe the integration points. 



appendix, we address the larger class of operators of twist two, of which the Konishi operator is 
the simplest (spin zero) representative. 

For the protected scalar operators (12. 2p the OPE takes the following general form: 

0{x,,y,)0{x,, y,) = T^CTX^T^^i^),, • • • {xu),s [o'^-'''''{x,) + ...]. (C.l) 

A,^,-^, (^12) ^ 

Here the sum on the right-hand side runs over conformal primary operators Q^'"^^'^ carrying 
Lorentz spin scaling dimension A and the dots denote the contribution of their conformal 
descendants. The relation ( IC.ll) generalizes ( 12. Sp which describes the most singular contribution 



of operators with the lowest value of A. Also, since each operator 0{xi,yi) belongs to the 
representation 20' of the R symmetry group SU{A), the right-hand side of (IC.ip involves the 
sum over all irreducible representations TZ that appear in the tensor product 

20' X 20' = 1 + 15 + 20' + 84 + 105 + 175 . (C.2) 

These representations can be identified from the ?/— dependence of the operators 0'^"''^^'^{x2) 
in the expansion. The contribution of each operator to the right-hand side of fIC.ip is accom- 
panied by the coefficient function Cqqq^. It determines the three-point correlation function 
{0{l)0{2)0'^"''^^'^{x3)) and depends, in general, on the coupling constant. 

A notable example of the operators - ^^igX shall play a special role in our discussion 
are the twist-two operator Os- They are SU{4) singlet bilinear operators with arbitrary (even) 
Lorentz spin S and naive scaling dimension A^°) = 2 + S. The Konishi operator is the special 
case of such operators with spin zero, S = 0. The scaling dimensions of the twist-two operators 
acquire an anomalous contribution: 



A5 = 5 + 2 + 75(a) = 5 + 2 + J]aSf , (C.3) 
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with 7^. being non-trivial functions of S. 

Applying (IC.ip . we find that every operator contributing to the OPE gives a definite contri- 
bution to the four-point function (12.31) known as a conformal partial wave amplitude (or CPWA) 



G(l, 2, 3, 4) = ^) ' (C.4) 

where G'^^\{u,v) describes the contribution of the conformal primary operator Q^'"^^'^ and 

its conformal descendants. The conformal partial waves G\ g{u,v) are definite functions of the 
conformal ratios u and v defined in f l2.13p . For m — )■ and f — )■ 1 they have the following 
asymptotic behaviour 



GaU% v) ~ - vf [1 + Oiu, l-v)] . (C.5) 

Expanding the known result for the four-point correlation function over the CPWA and applying 
(IC.Sp . we can extract the anomalous dimension of the operator. For more information see, for 
example [ill 115]. 

In A/" = 4 SYM, the twist-two operators are the ground states (or superconformal primaries) 
of "long" (or unprotected) supermultiplets. Each state (or superdescendant) in such a multiplet 
has different naive conformal dimension, Lorentz spin and SU (4) quantum numbers, but they all 
share the same anomalous dimension 7^(0). Therefore, to determine the anomalous dimension 
of the twist-two operators, we may look at the state in the multiplet which is most convenient 
to identify in the CPWA expansion ( ]C.4|) . As pointed out in |13], the best choice is the twist-six 



state corresponding to the S'f/(4) channel 105 in (1C.2|) . The advantage of this choice is that 
the twist-two supermultiplet has only one state in this SU (4) channel, while for all other choices 
there are multiple candidates^ 

The specific correlation function which singles out the state in the 105 has the form 

G,{1, 2, 3, 4) = {Oix,)Oix2)dixs)d{x,)) , (C.6) 

where O = ii{ZZ), O = tr{ZZ) and Z = $^ + is a complex scalar field. It can be obtained 
from the general expression for the correlation function fl2.3p by choosing the harmonic variables 
as Yi = Y2 = (l,i,0,0,0,0) and Y3 = Y4 = (1, -i, 0, 0, 0, 0). Then, the relations ([231) and (El 
take the following form 

^4(1,2,3,4) = ^^^^ X X ^ X 5^a^FW(x,) , (C.7) 

(47r2)4 xfsxii V ^ 

where the additional factor of u/v comes from the function -R(l, 2, 3, 4) (see Ref. [1]). Then, the 



^^Choosing the state in the 105 is helpful, but not indispensable for carrying out the CPWA analysis of the 
twist- two operators. The two- loop anomalous dimension of the S'C/(4) singlet state of spin two was found for the 
first time in [33], using a form of the CPWA expansion different from that in [33]. 
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CPWA expansion dCH) 

X ^F,{3 + S+^^sia),3 + S+^^sia),Q + 2S + ^sia)\v) . (C.8) 

where v = 1 — v and the sum runs over non-negative even spins S and the coefficient function 
^5(0) and the anomalous dimension 75(a) are functions of the 't Hooft couphng a 

Asia) = J2 , 75(a) = J2 <^'^f ■ (C.9) 

Comparing (1C.7P and ( 1C.8I) we thus obtain (for u ^ 0) 
HA ^'^^'^ = 5Z (1 - t;) t)^ 2F1 (3 + ^ + ^75, 3 + ^ + hs, 6 + 25 + 7s|^;) . 

e>l 5/26Z+ 

(C.IO) 

To make use of this relation, we expand both sides of this relation in powers of the coupling 
constant a and compare the coefficients of the powers of Inu. Furthermore, on the right-hand 
side, the contribution of the operator with spin S is suppressed by a factor of . Therefore, 
expanding both sides of (IC.lOP for w — )■ up to 0{v^) terms and equating the coefficients of the 
powers of Inw and v, we obtain equations for the expansion coefficients of As'{a) and 7s'(a) for 
S' < S. 

The main goal of this appendix is to apply (IC.lOp to extract the coefficients (1C.9P at three 
loops. To this end, we have to supplement (IC.lOp with the expressions for the correlation function 
p{£) fo^ £ = 1,2,3 in the OPE limit. At one and two loops, the corresponding expressions for 
F*^^) and F^'^\ Eq. (13. ip . are known analytically [30]. At three loops, we can use the result for 
F^^^ recently found in Ref. [1] (and described in section . The explicit expression for F^^^ in 
terms of the basis scalar integrals is 

=2(7(1,2,3,4) [xl,xl^h{l,2;3A) + xj^xj^ h{l , 3; 2 A) + x^xl, h{l A; 2 , 3)] 
+ 6[L(1,2;3,4) + L(l,3;2,4) + L(l,4;2,3)] 
+ 4[E(1;3,4;2) + F(l;2,4;3) + F(l;2,3;4)] 

+ (1 + 1/v) H{1, 2; 3, 4) + (1 + u/v) H{1, 3; 2, 4) + (1 + u) H{1, 4; 2, 3) . (C.ll) 

Here the integrals g, h were defined in (13. 2 p while the remaining integrals L, E and H can be 
found in [1]. The expressions for the integrals g, h, L are known in a closed analytic form. The 
integrals E and H on the other hand are not currently known analytically. However, for our 
purposes we only need their asymptotic behaviour at small u and the ffist few terms of their 
expansion in powers of v. 

The limit m — )■ 0, w — )■ can be described as the limit xi X2- For the E— and //—integrals, 
we apply the well-known formulae [55] for their asymptotic expansion in the limit Xi — )■ X2 typical 

^^Furthermore, by considering only the Hmit u — >■ with v finite, we are restricted to the sector where only the 
twist-two operators together with their superdescendants survive. Here it is known that there is a unique twist 
two, spin S supermultiplet for each even S. 
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of the correlation function flC.6P is (for m — )■ 



of Euclidean space written in terms of a sum over certain subgraphs of a given graplo- Making 
use of conformal invariance, we set 0:4 — )■ oo, xi = and arrive at the problem of analyzing the 
asymptotic behaviour of the Feynman integrals depending on the two external coordinates X2 
and X3 in the limit X2 — > 0. The conformal ratios take the form u = x\/x\ and v = x\.^/x\ so 
that for X2 — !■ px2 with p — > they scale in the Euclidean space as m = O(p^) and v = 0{p). We 
have evaluated terms up to order O(p^) for all integrals entering the right-hand side of (IC.lip 
and obtained the following results for the integrals 



+ 32^3 



I ( 1261 ^ _ 1060073 \ -A I ('_5_/- _ 1423 \ -3 _ ('_5_/- , - _ _3_/- 

~'~ V 576000 ^3 82944000/ ~'~V9216^3 110592/ V 1152^3 "T 3455; 128^3^ i6^3 

xj.xlMl, 4; 2, 3)=\nu [{§^^(3 + + (IC3 + ^)v' + (1^(3 + ^,)v' + i^Csv 

I ( 659 /■ _ 1682573 \ -,4 _ ( _65_^ , 1865 \ -,3 _ /J_a , .S^^J^S - _ 

~'~ V 144000^3 82944000/ I 9216^3 "T j^^^Qggg ; I 36^3 "T 3455 y 128 ^3 

2^132^24-^(15 2; 3, 4) = (2073600^^ + gife^^ + T^^^ + sli^ + ik) (^"^ 

+ ( 24679 ^4 _ ^2^^3 _ ^5^^2 _ ±^ _ A) V. 

~\ 4608000 18432 27648 64 64/^ "/ 

, ( 10888367 ^4 , ^0^^3 , i36L,;2 > ±y > ±) 

' \ 1244160000 ' 18432 ^ 82944 ^ 32 ^ 32 / 

I I 29383 A I 28271621 \„-A \ ( 5 a , 

V 16^5 "T 1728000^3 T 7464960000/ ~'~ V 64^5 ~'~ / 

+ (-iC5 + iic3 + MeV + (-iC5 + ^(3)^; ■ 



f _L iO]_/^ I 1381 \ -,3 
['^5 "T 4608*53 "T 331776 

1C5 + MC3-M' (C-12) 



and for the i7— integrals 



X 



X24-nyi-,C>, z,,^) — y 6912000*^ 110592*^ 10368*^ ^ 1536 ^ 192 M^^^ "/ 

, / 51643 -4 , 575 ^3 , ^1^-2 _ ±\ /j x2 

^ V 11520000 ^ 110592 ^ 41472 16/ V^^^"/ 

, ( 25959283 -4 _ 36013 -,3 _ 8623 -2 _ 15_- , _9_\ 1 

V 1244160000 1327104 248832 512 ~'~32/ 
I ( 1997 /■ I 380181271 \ -^4 /'__29_a , 344443 \; 

V 144000^3 "T 12441600000/ '\ 2304^3 "t" 7952624 ^ ' 

+ i-ieCs + ^y+{l2Cs + i2> + \Cs- 



2-132-24-" ll)^i^)<j; — U912000^ ^ 36864 ^ ^ 20736^ ^ 1536 ^^ 192;' "J 

, 1628267 ^4 _ 3155 ^3 _ ^7^^2 _ _ ±) \2 

^\ 69120000 110592 10368 64 16^ "/ 

, / 103801967 -4 , 136573 -,3 , 33173 -2 , _95_ , _9_\ 1 

~ V 1244160000 ~'~ 1327104''^ 248832 512 32/ " 

I ( 19253 A _ 282347651 \ -^4 , (119/- _ 928859 

V 144000 ^3 3110400000/ "'"^68 ^3 7962624 / 

_L /jra _ 119557 \ -,2 , (Ji^ _ 125\- , If _ J_ 
1432^3 746496^^ "'"^32^3 512/'^"'" 4^3 16 ' 

xl,xlH{l, 2; 3, 4) = Inn [(-^^(3 + jM^y + {-^3 + Ws)^' + (-^(3 + 7^)^' + U^] 

_L ( 217 /■ _ 791 ^-,4 I /J_A ^\ -3 I ^"1712 i A/- ,l _ 

V 11520^3 276480/ """^32^3 256/ """^128^3 1536 / "'"32^3'^ i6^3- 

(C.13) 



To find all the contributions to the asymptotic expansion in an automatic way we prefer to use the code 
asy.m |46) which reveals the contributions in the language of regions. Observe that this code works not only for 
limits typical of Euclidean space but also for other limits, in particular, for limits of Sudakov type. 
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Substituting these expressions into flC.10|) and going through the steps described above we 
obtain the following results for the three-loop anomalous dimensions 75(0) (for S = 0,2,4) 

21 

7o(a) = 3a - 3a^ + —a^ + O (a^) , 

, , 25 925 2 241325 3 ^ . 4. 

72 a = — a a H a + (J (a , 

' ^ ' 6 216 31104 ^ ^ ' 

, , 49 45619 2 300642097 3 ^ . 4^ , 

74(a) = —a H + O (an , (C.14) 

' 10 9000 32400000 ^ ' ^ ' 

and for the corresponding coefficients As{a) 

A,{a) = + + - (2C3 + ^ + 12) + O {a') , 

205 2/^5(3^ 76393 \ 3 / 1315(3 , 125(5 , 242613655 ^^ , ^ ^ 4^ 

7^2 a = —a h a — a + (a , 

^ ^ 1764 V 28 148176 J \ 5292 168 112021056 J ^ ^ ' 

, , , 553 2 7(3 880821373 \ 
AAa) = —a— h a h 



54450 V 440 17249760000 / 

3 



a 



/520093C3^35C, 1364275757197X 

V 26136000 528 5692420800000/ ^ ^ ^ ^ 



We verified that the relations (1C.14P are in agreement with the results of [24J. Notice that 



7o(ci) coincides with the anomalous dimension of the Konishi operator 7x;('3.), Eq. (II. ip . and the 
expansion coefficients of Ao(ci) = '^i^i'^^^^o,^ coincide with the coefficients a'^^^ in (I2.16p . 

Finally, we should note that it is possible to invert this entire process, namely to use predic- 
tions for the normalization and anomalous dimensions of operators, and plug them into the con- 
formal partial wave expansion (IC.4P to obtain predictions for the four-point correlation function. 
More specifcally, by using the predicted three-loop, all-spin, twist-two anomalous dimensions 
of [24], together with two- loop and one- loop twist-two data, and plugging into ( ]C.8p . one can 
obtain the three-loop correlation function in the limit m — )■ with v finite (all except the finite 
part as M — )■ which would require three- loop normalization). Summing the resulting expansion, 
one obtains a closed analytic form for the correlation function in this limit as a sum of harmonic 
polylogarithms with argument v 
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